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ABSTRACT 


The  fundamental  concept  of  orthogonality  of  mathematical  objects  occurs  in  a  wide 
variety  of  physical  and  engineering  disciplines.  The  theory  of  orthogonal  functions,  for 
example,  is  central  to  the  development  of  Fourier  series  and  wavelets,  essential  for  signal 
processing.  In  particular,  various  families  of  classical  orthogonal  polynomials  have  tra> 
ditionally  been  applied  to  fields  such  as  electrostatics,  numerical  analysis,  and  many 
others. 

This  thesis  develops  the  main  ideas  necessary  for  understanding  the  classical  theory 
of  orthogonal  polynomials.  Special  emphasis  is  given  to  the  Jacobi  polynomials  and  to 
certain  important  subclasses  and  generalizations,  some  recently  discovered.  Using  the 
theory  of  hypergeometric  power  series  and  their  q  —extensions,  various  structural  prop¬ 
erties  and  relations  between  these  classes  are  systematically  investigated.  Recently,  these 
classes  have  found  significant  applications  in  coding  theory  and  the  study  of  angular 
momentum,  and  hold  much  promise  for  future  applications. 
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I.  INTRODUCTION 


The  abstract  concept  of  orthogonality  of  functions  (or  other  mathematical  objects) 
is  a  generalization  of  the  notion  of  having  two  or  more  vectors  perpendicular  to  one 
another.  This  concept  arises  naturally  in  a  wide  variety  of  physical  and  engineering 
disciplines.  For  example,  the  theory  of  ortliogonal  Ametions  is  central  to  the  develop* 
ment  of  Fourier  series  and  wavelets  which  are  essential  to  signal  processing. 

Classical  Fourier  series  (real  form)  depend  on  the  property  that  the  trigonometric 
Ametions  sine  and  cosine  arc  orthogonal  (on  an  appropriate  real  interval)  in  a  formal 
sense  that  will  be  made  precise  later.  As  a  consequence,  a  bounded  periodic  Ainction 
f{x)  of  period  2ir  which  satisfies  the  Dirichlet  conditions!  may  be  expressed  in  the  fonn 

OO 

fix)  y  Oo  +  ^  cos  /IX  +  6,  sin  /tx ) 


where 


I  /W  cos/ixiix,  /I  •0,1,2,... 

I  /W  sin/fxrfx,  /I- 1,2,3,... 

are  the  classical  Fourier  co^cients.  These  formulas  can  be  modiAcd  via  a  change  of 

variable  to  accomodate  any  such  function  of  period  2L.  [Ref.  1;  p.  529) 

This  property  can  be  used  to  generate  other  classes  of  orthogonal  functions  • 

polynomials,  for  example  -  that  behave  in  very  structured  and  useful  ways  such  as  in 

generalized  Fourier  series.  In  particular,  specific  families  of  these  "classical  orthogonal 

polynomials"  have  traditionally  been  used  for  solving  problems  arising  in  various  areas 

* 

of  applied  mathematics,  physics,  and  engineering,  among  others. 

This  thesis  develops  the  main  ideas  necessary  for  understanding  the  classical  theor}' 
of  orthogonal  polynomials.  Special  emphasis  is  given  to  the  Jacobi  polynomials  and  to 
certain  important  subclasses  and  generalizations.  Much  of  the  investigation  will  be 

1  Dirchlut  conditions:  (i)  In  any  period /(x)  is  continuous,  e.xccpt  possibly  for  a  finite  number 
of  jump  discontinuities,  (ii)  In  any  period /(x)  has  only  a  finite  number  of  maxima  and  ininitnu. 
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made  using  tlie  theory  of  hypergeometric  power  series  and  their  q  -extensions.  The 
classes  discussed  in  this  thesis  arc  but  a  small  fraction  of  those  identified  and  studied  in 
the  literature. 

A.  CHEBYSHEV  POLYNOMIALS 

The  Chtbyshev  polynomials  qf  the  first  him/,  TJix),  arise  from  an  elementary  trigono* 
metric  consideration.  As  such,  they  satisfy  various  properties  and  identities  which  are 
easily  derived  directly  from  their  definition,  many  of  which  are  observable  from  their 
graphs  (see  below).  This  class  of  polynomials  will  serve  as  the  model  for  some  of  the 
basic  smicture  of  more  general  classes. 

For  /f- 0,1,2, ... ,  define 

r„(.t)  -  cos(rt  arccos  x),  —  1  SI  x  S  1 

i.e.,  letting x ■> ccs 0 ,  0^6 ^n, 

(1)  r„(  cos  9)  ■>  cos  ftO. 

Some  immediate  consequences  of  (1)  are  |r,(x)|^l  for  |x|  1,  with 

cos ^  “ ( -1)*,  0^k^n  \  in  particular  r,(l) »  1  and  Tj[  -1) - ( -1)"  for  all  n 
which  can  be  seen  graphically  in  Figure  1. 

1.  Three'term  Recurrence  Relation/DlfTerential  Equation 
From  (1),  we  have 

(2)  7'o(x)-l  and  r,(x)-x, 
and  by  considering  the  identity 

(.1)  cos(a  +  i)  +  cos(a  -  fc)  •  2  cos  a  cos  b 

with  a"*n0,b^0,  we  obtain 

(•<)  W  ■  2jf7;(x)  -  T^iix). 

Equation  (4)  is  known  as  the  three-term  recurrence  relation  for  TJ[x)  which  to¬ 
gether  with  initial  conditions  (2)  imply  that  r,(x)  is  a  polynomial  of  degree  exactly  n, 
called  the  n*  Chebyshev  polynomial  of  the  first  kind.  Note  that  the  leading  cocllicicnt 
of  r,(x)  is  2^'  for  n'^2.  An  inductive  argument  applied  to  this  recursion  shows  that 
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Figure  1.  Chebyshev  Polynomials 


TJipc)  is  an  even  function  if  a  is  even,  and  odd  if  a  is  odd  (see  Figure  1).  The  first  few  are 
listed  below: 

1  •  7',(x)  -  X ,  r2(x)  -  2x*  -  I ,  TjCx)  -  4x^  -  3x 

r4(x)  -  8x*  -  8x*  +  I ,  Tsix)  -  I6x*  -  20x^  +  5x. 

Diiferentiating  (1)  twice  with  respect  to  B  yields  the  second  order  differential 
equation  for  r,(x) : 

(3)  (1  -  -  xTM  +  -  o. 


3 


4 


(7) 


1 


i:- 


WT,(x) 


'dx  —  h„^  8 


m,rt* 


This  important  property  is  formally  known  as  the  orthogonality  relation  for  the 
Chebyshev  polynomials.  The  reason  for  this  temunology  will  become  clear  in  the  next 
chapter. 

3.  Zeros  of  Chebyshev  Polynomials 

2k  —  \ 

Setting  T„(  cos  $)  =»  cos  nd  =  0,  we  obtain  0  =  0*,,  =  — ^ —  n,  ,  i.e., 

X  =  jr*, ,  =  cos  0*. , ,  1  ^  ^ 

Thus  all  the  zeros  of  T„(x)  are  real,  distinct,  and  may  be  regarded  as  the 
projections  onto  the  interval  (—1,1)  of  the  equally  distributed  points  0*,,  on  the  unit 
circle,  as  seen  in  Figure  2.  Moreover,  for  1  ^  A:  ^  an  easy  algebraic  check  verifies  that 
<  ^*+1. and  therefore  <  Jf*.,  <  Hence,  the  zeros  of  T^,(x) 

interlace  with  those  of  T„(x).  This  interlacing  of  zeros  is  a  striking  feature  of  the  plots 
in  Figure  1. 

The  zeros  of  Chebyshev  polynomials,  and  of  other  orthogonal  polynomials  in 
general,  are  extremely  important  for  applications  to  numerical  analysis,  electrostatics, 
and  many  other  fields. 

4.  Looking  Ahead 

Many  of  the  properties  derived  for  the  Chebyshev  polynomials  T„(x)  from  their 
trigonometric  definition  (1),  extend  to  more  general  classes  of  orthogonal  polynomials 
via  a  general  theory,  elements  of  which  will  be  devciOped  in  this  thesis.  Some  of  the 
many  properties  satisfied  by  these  classes  that  we  will  derive  include: 

1.  Orthogonality  with  respect  to  a  weight  function 

2.  Three-term  recurrence  relation 

3.  Second  order  differential  or  difference  equation 

4.  Hypergeometric  series  expression 

5.  Rodrigues'  formula 

6.  Generating  function. 

The  general  approach  we  will  take  is  to  define  these  "classical  orthogonal  polynomials" 
via  terminating  hypergcometric  power  series,  and  from  this  prove  (most  ol)  the  other 
properties.  However,  because  of  this  equivalence,  many  authors  choose  to  define  a  given 
class  using  one  of  these  other  characterizing  properties. 
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Figure  2.  Zeros  of  Chebyshev  Polynomials 


In  order  to  understand  the  general  theory,  it  is  first  necessary  to  define  the  ab¬ 
stract  concept  of  orthogonality  in  an  appropriately  defined  "space*  of  functions.  We 
turn  our  attention  to  these  fundamental  ideas  in  the  next  chapter. 
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II.  BACKGROUND 


A.  ELEMENTARY  LINEAR  ALGEBRA 

1.  Vector  Spaces 

Let  R"  denote  the  collection  of  all  vectors  (n-tuples),  UB((i|,a2, ...  ,0;,),  where 
each  a,  e  R,  /  >■  1, 2, ... ,  n.  The  standard  Euclidean  inner  product  (also  referred  to  as  the 
dot  product )  of  two  such  vectors  u  =  (<*1,  Oj, ... ,  <^,)  and  v  =  (6„  bj, ... ,  b„)  is  given  by 

n 

<u,v> 

iml 

The  length,  or  norm,  of  a  vector  u  6  R"  is  given  by 

Hull  =  V<u,u>  =  I  a,|  ^ 

Two  vectors  u,  v  e  R"  are  perpendicular,  or  orthogonal,  if  and  only  if  <  u,  v>  =0. 

The  objective  of  this  chapter  is  to  extend  these  familiar  notions  to  objects  other 
than  classical  Euclidean  vectors,  in  particular,  the  “vector  space"  of  polynomials  defined 
on  a  real  interval  [a, 

A  vector  space  V  over  a  scalar  field  F  (usually  R  or  C  )  is  a  nonempty  set  of 
objects  called  vectors,  for  which  the  operations  of  addition  and  scalar  multiplication  are 
defined.  Addition  is  a  rule  for  associating  with  each  pair  of  vectors  u  and  v  in  an 
element  u  +  v,  called  the  sum  of  u  and  v  .  Scalar  multiplication  Is  a  rule  for  associating 
with  each  scalar  c  in  f  and  each  vector  u  in  K  an  element  cu,  called  the  scalar  multiple 
of  u  by  c.  [Ref.  2:  p.  150] 

For  all  u,  V,  w  6  V  and  c,de  F,  a.  vector  space  V  must  satisfy: 

1.  Additive  closure,  u, ve  Fa^u  +  ve  V 

2.  Commutativity,  m-  v  *  v  +  u 

3.  Associativity,  u  +  (v  +  w)  =  (u  +  v)  +  w 

4.  Additive  identity.  There  exists  a  zero  vector,  Oe  V,  such  that  0  +  u  =  u  4-  0  =  u. 

5.  Additive  inverse.  For  each  u  6  V,  there  exists  a  vector  -  u  e  F,  such  that 
u  4-  ( -  u)  =»  ( -  ii)  +  u  =  0. 

6.  .Multiplicative  closure,  u  e  F  and  ce  F^  cue  V 
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7.  Distributivity.  c(u  4- v)  —  cu  +  cv 

8.  Distributivity.  (c  +  <^0  "  cu  + 

9.  Multiplicative  associativity.  c(</u)  *  (c<Ou 

10.  Multiplicative  identity.  There  exists  a  scalar  1  e  F  such  that  lu  =  u. 

Example  1:  R”  (the  model) 

Example  2:  «  (polynomials  of  degree  on  the  interval  Ca,i>]) 

Example  3:  PCa,6]  -  (polynomials  on  the  interval  [a, 6]} 

Example  4:  Cta,6]  -  (continuous  functions  on  the  interval  [a, 6]} 

Note  that  PuiOjbl  c  P[a,b1  a  C[a,6].  These  vector  spaces  are  sometimes  re¬ 
ferred  to  as  function  spaces.  The  interval  may  be  finite  or  infinite  (i.  e., 
[a,  oo) ,  (— oo,  A] ,  or  (-00,  oo) )  for  our  purposes. 

A  subset  1/  of  a  vector  space  V  is  said  to  be  a  vector  subspace  of  V  if  it  is  a 
vector  space  in  its  own  right. 

Example  5:  is  a  subspace  of  PCa,6],  which  in  turn  is  a  subspace  of 

C[a,£>]. 

Given  a  set  of  vectors  (v„  Vj, ... ,  v,}  in  a  vector  space  V,  and  scalars 

Cl,  c, . .  the  vector  c,Vi  4-  CjV,  H - V  c,v„  e  is  said  to  be  a  linear  combination  of 

(v„  Vj, ... ,  vj.  The  set  of  vectors  (v„  Vj, ... ,  v,}  is  said  to  be  linearly  dependent  if  there 
exist  scalars  c„  Cj, ... ,  c„  not  all  equal  to  zero,  such  that  the  linear  combination 
c,v,  4-  CjV,  -f-  •••  4-  c,v„  =  0.  (Equivalently,  at  least  one  of  the  vectors  v,  can  be  expressco 
as  a  linear  combination  of  the  others.)  Otherwise,  the  set  (v„  v, . v,}  is  linearly  inde¬ 

pendent.  An  infinite  set  S  —  (v„  v„ ... ,  v„ ... }  is  defined  to  be  linearly  independent  if  every 
finite  subset  of  S  is  linearly  independent;  otherwise  S  is  linearly  dependent  (Ref.  3:  p.  8]. 
The  vectors  (y„  v„ ... ,  vJ  are  said  to  span  V  if  every  vector  v  e  K  can  be  represented  as 
a  linear  combination  of  (v„  v„ ... ,  vJ.  In  this  case,  we  write  V  ■»  span{\x,  v„ ... ,  v„}.  The 
vectors  (v„  Vj, ... ,  >,}  form  a  basis  for  V  if  they  are  linearly  independent  and  span  V.  The 
dimension  of  K  is  the  number  of  elements  in  any  basis. 

Example  6:  The  set  (e,,  e„ ... ,  e,}  is  the  standard  basis  for  R",  where 
e,  =»(0, 0 . 0, 1,0,...  ,0)  i.e.,  the  vector  with  a  one  in  the  /'*  position  and  zeros  else¬ 

where,  1, 2, ...  ,n. 

Example  7:  The  set  (1,  jr,  .x^...  is  the  .standard  basis  for  PJ.a.b'].  (Linear 
independence  is  ensured  by  the  Fundamental  Theorem  of  Algebra.)  The  dimension  of 
Pi,lo,b']  is  therefore  A4-1. 
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Example  S:  The  set  ...  ,jc*, ... }  is  the  standard  basis  for  P[a,b'\:  and 

hence  P[a,6]  is  an  infinite-dimensional  vector  space. 

2.  Inner  Product  Spaces 

An  inner  product  on  a  real  vector  space  K  is  a  mapping 

<  ,  >:rx  K-»R 

such  that  for  all  u,  v,  w  e  K  and  oe,  /?  e  R,  the  following  properties  hold: 

1.  Positive  definiteness:  <  u,  u>  2  0,  and  <  u,  u>  0  if  and  only  if  u  «  0 

2.  Symmetry:  <u,y>  —  <v,  u> 

3.  Bilinearity:  <att  +  j8v,  w>  ■■«<u,w>  +)J<v,w> 

A  vector  space  with  an  inner  product  is  known  as  an  inner  product  space. 

Example  9:  R* ;  let  constant  'weights'*  w,  >  0  be  given,  /  —  1, 2, ... ,  n. 

For  u  -  (a„  . and  v  -  (h,.  6,. ... ,  6,),  u,  v  e  F, 

n 

<u,v> 

im\ 

If  w,*  1  for  I'o  1, 2, ...  ,ff,  then  this  reduces  to  the  standard  Euclidean  inner  product, 
or  dot  product.  Otherwise,  this  is  referred  to  as  a  weighted  inner  product. 

'The  next  two  examples  are  commonly  applied  inner  products  on  function  space, 
and  are  analogues  of  the  previous  example.  We  assume  a  given  weight  function 
w(jf)  >  0  in  (a,6),  integrable  in  the  first  case  (c.g.,  continuous  for  [o,f>]  a  finite  interval). 
Example  10:  K-  PJia,b'\,  ?[«,/>],  C[«,63 


<f,g> 


Example  11: 


N 

-  2]/Wg(x)w(jf) 
x-0 

(Positive  definiteness  is  ensured  by  the  Fundamental  Theorem  of  Algebra.) 
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The  norm  induced  by  the  inner  product  is  given  by  ||u||  -  J<u,u>  .2 
Example  12:  For  the  inner  products  of  Examples  10  and  1 !  therefore 

ll/ll -(|^/(x)  ]’»<*)<&)''’  and 

respectively.  These  are  sometimes  referred  to  as  'Z,*-norms.* 

Two  vectors  u,  v  e  F  are  said  to  be  orthogonal,  denoted  ulv,  if  and  only  if 
<u,v>  —0.  The  vectors  u  and  v  are  said  to  be  orthonormal  if  ui.v  and 
||u|| » ||v||  >■  1.  Note  that  the  orthogonality  of  vectors  in  a  space  is  determined  by  the  in¬ 
ner  product  being  used. 

The  two  e.xamples  which  follow  refer  back  to  Chapter  1,  Section  A.l. 

Example  13:  Formula  (6)  shows  that  the  functions  { 1,  cosx,  cos  2x, ... }  are 
orthogonal  on  [0,  Jt]  with  respect  to  the  uniform  weight  function  w(x)«»  1.  A  similar 
computation  shows  that  the  same  property  holds  on  [-»r,  jt]  with  respect  to  the  weight 
function  m<x)  « ,  i.e., 


^  >  =4“  f  fix)g{x)dx. 


One  advantage  of  preferring  this  inner  product  over  the  standard  one  lies  in  the  com¬ 
putation  of  norms.  Using  ||/||  ,  we  have  ||  1 1|  -  2  and  ||  cos  nx  ||  -  1  if 

n  S  1.  Hence  the  functions  {  1/2,  cosx,  cos  2x, ... }  are  orthonormal  on  [-«,  w]  with 
respect  to  the  inner  product  above.  Similar  statements  hold  for  the  integral  of  a  product 
of  two  sine  functions  on  [-«,  r],  as  well  as  for  the  product  of  a  sine  and  a  cosine. 

Example  14:  By  (7)  the  Chebyshev  polynomials  {7i,(x)}  form  an  orthogonal 
class  with  respect  to  the  inner  product  of  Example  10  above  on  [  -1,1]  with  the  weight 
function  w{x)  =  (1  —  x’)'*'*. 

2  Recall  that  <  u,  u  >  ^  0.  We  remark  that  in  the  same  way  we  defined  inner  product  earlier, 
it  is  possible  to  dcfuie  a  general  norm  on  a  vector  space  wliich  is  not  induced  by  an  iimcr  product. 
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We  remark  here  that  Examples  10  and  11  can  be  unified  into  a  single  inner 
product  on  a  "polynomiar  space  V  via 


fix)  gix)  rfa(x), 

where  </a(jr)  is  a  positive  Lebesgue-Stieltjes  measure  on  a  measurable  set  £  possessing 

finite  moments,  i.eM  jr  d’a(jr)  integrable,  a  ■>  0. 1, 2 . In  Example  10,  £  ~  [a, 6]  c  R  and 

doi{x) « iv(jr)<6:;  the  resulting  expression  is  known  as  a  continuous  inner  product,  while 
in  Example  11  the  set  £  consists  of  a  finite  number  of  points  (0, 1, ... ,  /Y}  cR,  and  the 
associated  .measure  gives  rise  to  a  discrete  inner  product. 

B.  FOURIER  SERIES 

Let  yeV,  and  (/  be  an  /t-dimensional  subspace  of  V  having  some  orthonormal  basis 
( u„ ... ,  u„ }.  (Any  basis  can  be  orthonormalized  via  the  Gram-Schmidt  process  -  see  next 
section.)  The  vector  v  can  be  resolved  into  a  sum  of  two  vector  components: 

(1)  v  -  (v  —  w)  +  w 

where  we  U  and  (v-w)l  £/,  (See  Figuio  3.)  The  vector  w  is  referred  to  as  the 
orthogonal  projection  of  v  onto  (/.  Since  the  vector  (v  —  w)  is  orthogonal  to  every  vector 
in  U  by  construction,  it  follows  that  for  each y  —  1, 2, ... , »,  <  v-w ,  u,  >  »  0,  or 

(2)  <  V,  U;  >  ■»  <  W,  Uy>. 

Moreover,  since  it  lies  in  U,  vector  w  can  be  expressed  as  some  linear  combination  of 

{ u„  ... ,  n, } : 

n 

w-]^c,n,. 

i-i 

Take  the  inner  product  of  both  sides  with  Uy  for  eachy  >*1,2,...,  n.  From  the  assumption 
that  <  u, ,  Uy  >  -0  unless  i  -y,  we  have  the  property  that 


via  (2)  and  the  assumption  that  <  >  >  |  f  *  1. 

Thus, 


(4)  w-2<v,u,>u, 

and  this  vector  represents  the  '1)est  approxinuition"  in  C/  to  v  «  K'  in  the  sense  that  of 
all  vectors  z  e  (/,  it  is  the  projection  vector  w  e  {/  which  uniquely  minimizes  the  distance 
Bt-xll. 

% 

Suppose  now  that  U  is  an  inrinite<dimensional  subspace  of  V  (also  infinite* 

•  dimensional),  having  orthonormal  basis  { U|, ... ,  it,, ... }.  Then  from  (1)  and  (4) 

A 

T-(V-W)  +  ^  <  V,  U,>  U, 

1-1 

we  may  write 

00 

(5) 

1-1 

in  the  sense  that 

(6)  lim  II  v-w  II  -  0, 

A  '•OO 

i.e.,  the  norm  of  the  "residual  vector"  (and  hence  the  vector  itseU)  v-w  -+  0  as 
a  “*  OO.  Formula  (5)  is  known  as  the  generalatd  Fourier  series  for  v  e  K  with  respect 
to  the  orthonormal  basis  { u, The  coefllcients  given  in  (3)  are  called  the  fe/»cra/«e</ 
Fourier  coefficients  of  v  e  V.  Statement  (6)  is  known  as  the  norm  convergence  property 
of  Fourier  series,  and  the  "minimization  property"  mentioned  above  extends  to  this 
infinite-dimensional  case. 

Example  15:  Let  F— Cta,6],  and  be  an  orthonormal  basis  of 

eigenfunctions  (sometimes  refen’cd  to  simply  as  an  eigenbasis  )  of  V.  Then  /  e  V  has  a 
Fourier  series  representation 


OO 


1-0 


12 


Flgurt  3.  B«st  Approximation 


with  Fourier  coeflicients 


(7)  c,  -  <  -  f  fix)  ^Xx)  Mx)  dx. 

In  this  Ainction  space  context,  norm  convergence 


Urn 


0 
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is  reierred  to  as  mean  square  convergence,  and  is  the  least  squares  principle  in  regression 
analysis. 

In  particular,  if  V  is  equipped  with  the  inner  product  of  Example  13  and 
orthonormal  basis 

{ }  “  { 1/2,  cos  X,  cos  2x, ... ,  sin  x,  sin  2x, ... } 

on  [-R, «]  (see  Example  13,  Section  A.2),  then  a  suitable  flmction  /•  V  (and  its 
2>t«periodic  extension  on  R  )  has  a  classical  Fourier  series 


OO 

/(jf) y 2 cosnx  +  A„  sin«x) 

naO 


where 


</,  cos/uc>  j  f{x)  cosnxdx 
</»  sin«jr>  J  f{x)  svanxdx, 

as  indicated  in  the  Introduction. 

C  GRAM-SCHMIDT  ORTHONORMALIZATION 

The  Gram>Schmidt  process  orthonormalizes  any  set  of  linearly  independent  vectors 
in  an  inn<  r  produc,  space.  This  method  will  be  used  in  later  sections  for  dilTercnt  inner 
products  on  tne  vector  space  P[a,6j. 

Begin  with  an  inner  product  space  V  and  any  set  of  vectors  {v„  v„ ... ,  v„ ... ),  finite 
or  infinite,  such  that  any  finite  number  of  elements  of  this  set  are  linearly  independent. 
Recursively  define  a  new  set  of  \  ectors  {u„  u„ ... ,  u,, ... } 


lly*ll  ’ 


*-1.2, 


a,,.. 


where  y*  -  v*  -  w*,  with3 


0 

3  By  convention,  *  0,  giving  Wi  =«  0. 

Cl 
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These  new  vectors  {u„  Uj, ... ,  u., ... }  are  orthonormal  by  construction  and  span  the 
same  space  as  the  original  vecton.  Note  that  this  process  occurs  in  two  stages: 
orthogonalization  and  normalization.  The  orthogonalization  is  accomplished  by  sub¬ 
tracting  w«,  the  orthogonal  projection  of  v«  onto  the  subspace  spanned  by 
{U|,  U), ... ,  The  component  of  which  remains,  denoted  above  as  yj^,  is  tlicn 
orthogonal  to  the  vectors  (u,.  U], ... ,  u^,)  as  shown  in  Figure  4.  The  normalization  is 
then  achieved  by  dividing  y^  by  its  norm,  thus  giving  it  unit  "length*’. 

1.  Legendre  Polynomials 

Example  16:  Let  I'"  PC -1,1]  with  basis  {l,jf,x*, ...  ,ac*, ... }  and  uniform 
weight  function  w{x)  «■  1.  The  inner  product  is  then  given  by  <f,g>  “  f\f{x)gix)<dx. 
The  Gram-Schmidt  process  yields  the  set 

as  an  orthonormal  basis  for  PC-1,1].  Since  this  set  is  linearly  independent,  we  can 
standardize  the  set  by  taking  scalar  multiples  of  these  polynomials  so  that  P,(l)»  1. 
Members  of  the  resulting  orthogonal  set 

( W)Z.  -  |l.  *.  y  (3x>  -  1),  i  (ir*  -  ix), ... } 

are  known  as  the  Legendre  polynomials  on  [-1,1].  If  the  normalized  Legendre 
polynomials  { p^x)  }i!!«  are  used  as  tlte  orthogonal  eigenbasis  for  a  Fourier  series,  the  re¬ 
sulting  expansion  is  oflen  referred  to  as  a  Legendre  series  reprc.sentation;  when 
Chebyshev  polynomials  are  used,  we  obtain  a  Chebyshev  series  representation,  etc. 

The  Gram-Schmidt  process  can  always  be  used  in  this  way  to  generate  a  class 
of  orthogonal  polynomials  with  respect  to  a  given  inner  product  (i.e.,  weight  function) 
on  a  real  interval.  When  using  the  Gram-Schmidt  process  from  the  basis 
{l,x,  jr*, ...  ,x", ... },  the  orthogonalization  stage  producing  y^,  results  in  a  .set  ol'  monk 
polynomials,  i.e.,  die  leading  cocdlcicnt  of  each  polynomial  is  one.  In  the  normalization 
stage,  we  are  dividing  by  the  norm  ||  y«  ||  >  0.  Thus  the  loading  coefficient  of  polynomials 
in  an  orthogonal  class  is  strictly  positive.  In  die  next  chapter,  wo  will  examine  other 
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ways  to  define  these  classes.  It  is  the  structure  and  applications  oF  certain  of  these 
classes  with  udiich  we  will  primarily  be  concerned. 

D.  THE  GAMMA  FUNCTION 

The  gamma  Jmction  r(x)  is  a  flmdamental  mathematical  object  that  appears  fre* 
quently  in  the  represenutions  of  orthogonal  polynomials  as  well  as  in  many  other  ap¬ 
plications.  This  'special  fimction*  was  developed  as  a  generalization  of  the  factorial 
function  of  the  natural  numbers.  As  we  will  see,  the  gairuna  function  has  the  value 
(/f  -  1)!  for  the  positive  integers  n  but  it  is  defined  for  noninteger  values  as  well. 

A  conventional  definition  for  the  gamma  function  is 
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The  positivity  ofx  ensures  that  this  improper  integral  converges.  We  now  develop  some 
fundamental  properties  of  the  gamma  function.  Integration  by  parts  in  (S)  yields 

(9)  r(jf+I)-xr(x). 

We  now  introduce  the  Poehhammtt  ^mbot  or  factorial,  (a), ,  to  simplify  our 
notation.  For  it  >  0.  define 

(a)„  ■  <i(a  +  IXa  +  2) ...  (a  +  It  -  1)  ,  if  a  >  1 

and  (a)«  >■  1.  Letting  a  -•  1  gives  (1),  ■  (1)(2)(3) ...  (a)  -  a!.  Note  that  for  a  negative  inte¬ 
ger,  (— ai),  —  0  if  a>at>0.  The  shifted  factorial  can  be  defined  for  negative  subscripts 
but  we  will  not  need  this  in  our  work  with  polynomials. 

Iteration  of  (9)  a  times  yields 


(10)  r(x  +  a)-(x),r(x) 

for  every  positive  integer  a.  Using  this  property,  the  gamma  function  can  be  extended 
to  include  negative  real  numbers  by  defining 

(11)  r(jf)  —  "pr— r(j:  +  a)  for  — a<x<“a-l-l. 

Since  this  expression  is  undefmed  when  x  is  zero  or  a  negative  integer,  the  gamma 
function  is  not  defined  for  those  values. 

Letting  x«  1  in  (8)  and  computing  directly,  we  have  r(l) 1.  It  then  follows  that 
r(a  -I- 1)  •  a!  by  letting  x  «  1  in  ( 10).  Furthermore, 

^0  Jo 

where  the  second  integral  can  be  evaluated  by  standard  methods  involving  multiple  in¬ 
tegrals. 


Finatty,  we  define  a  generalized  binomial  coejfieient  as  follows.  For  jr  and  «  non* 
negative  integers,  define 


For  nonintegral «,  define 


.xT  (•);,  "  r(x+i)r(«+i)  • 


1.  The  Beta  Function 

An  integral  related  to  the  gamma  function  defines  another  usefbl  function  called 
the  beta  function  which  is  given  by 


(12) 


B{x,y) 


i 


for  x,;^>0.  We  now  establish  an  important  connection  between  the  beta  and  gamma 
functions.  We  start  with  an  identity  easily  verified  from  (12): 

(13)  ^Jr,;.+  l)-^x,jr)-^x+l,j.). 


Also  from  (12) 


B(x+l.>r)- 

which  when  integrated  by  parts  gives 

B(x+  1,  y)^yBix,  1). 
Substituting  into  (13),  we  obtain 

Bix,  J/+  1) 

which  when  iterated  yields 


IS 


.  (Jf+A  .  (Jf+A 


f  (’-'(i-(r'“'<ft. 

vn 


Changing  variables  from  /  to  ~ , 


(14)  s(^,j,).iL^fV-(i-i-)«'-V,. 

(A«  Jo  V  "  / 

Taking  the  limit  as  n  -*  cxs,  and  using  the  fact  that  lin}(l-  we  have 

(15)  »(*.5')  =  n^)l!S,7T^- 

(A  n 

(The  fact  that  we  can  pass  the  limit  through  the  integral  on  the  right  can  be  mathemat¬ 
ically  justified.)  If  **  1 ,  then  (15)  gives 

B(x,I)-r(Ar)ltaiif^. 

'»-*«>  nl  n 


By  direct  evaluation  using  (12), 


c,l)=  I  = 

Jf\ 


Hence 


-V  'n-*oo  n\n* 


which  can  be  written  as 


r(.Y)  =  lim-7^-f-:7-. 

rt-oo  jr(;f  +  1)„ 


Noting  that  .x(Ar  4- 1)„  =  (.y),+,  =  (.Y)„(jf  +  n),  we  have 
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which  gives  the  form 


«-.oo  (x)„(x+n) 


(16)  r(jc) » lim  . 

'  n-oo  (jf)„ 

(Equation  (16)  was  Euler's  original  definition  of  the  gamma  function.  A  separate  "esti¬ 
mation'  argument  may  be  used  to  show  that  this  limit  mathematically  exists.)  Thus  by 
(15). 

ix+y)n 


Then  by  (16),  we  have  the  useful  identity 


(17) 


Bix,  y) 


rjxWiy) 

r{x+y)  • 


We  will  find  this  identity  useful  in  understanding  the  Jacobi  polynomials  in 
Chapter  IV,  where  it  becomes  necessary  to  evaluate  a  related  integral: 


f’  (l-jc)*(l+x)^</x. 

We  remark  here  for  future  reference  that  the  formal  change  of  variable  x  *=  1-2/  can  be 
used  to  transform  this  integral  into 


2«+/i+i  J '  j ^  2*+^+'  5(a+ 1,  p+l) 

^  r(«4-i)r(/?-n) 

ria+P+2) 
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III.  GENERAL  THEORY  OF  CLASSICAL  ORTHOGONAL 

POLYNOMIALS 

In  this  chapter,  we  examine  some  of  the  characteristic  properties  associated  with 
classes  of  orthogonal  polynomials.  Some  of  these  properties  provide  alternate  means 
of  deflning  a  class.  These  alternate  definitions  often  provide  a  straightforward  way  of 
producing  a  specific  result  that  may  be  very  difficult  to  derive  othen^ise. 

Throughout  this  chapter,  we  let  {/>„(x)}“o  denote  a  set  of  real  polynomials  with 
p„{x)  of  degree  n.  i.e.,. 

(1)  Pnix)^Kx''  +  +•••  ,  A„>0. 

Recall  that  these  polynomials  are  said  to  be  orthogonal  on  an  interval  [a,f>]  with  respect 
to  a  continuous  weight  function  w(jf)  >  0  on  {a,b)  if 

(2)  <  Pm^Pn>  Pmi.x)Pn{x)Mx)dx  -  /l"' 

•<1 

where  the  normalization  A‘‘  vt  0  is  chosen  to  simplify  the  expression  of  certain  formulas. 
Note  that  since 

II  Pn  11^  -  f  [  Pnix)'fwl.x)<ix  ~h~\ 


it  follows  that  A  >  0. 

A.  POLYNOMIAL  EXPANSIONS  . 

We  begin  by  showing  that  any  real  polynomial  qjix)  of  degree  m  on  [a,f>]  can  be 
written  as  a  linear  combination  of  orthogonal  polynomials  { p,(jf)}7.o : 

m 

(3)  = 

laU 

for  constants  tx,^ ,  /  =  0,  1,  2, ... ,  m. 

The  proof  is  by  induction  on  the  degree  m.  Since  qj{x)  is  a  real  polynomial,  wc 
write 
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where  a„  ^  0.  For  m  -  0,  (3)  reduces  to  fn,  -  ao.o  using  the  form  for  po{x)  given  in  (1) 
and  ffg,,  is  uniquely  determined. 

For  the  induction  hypothesis,  suppose  tliat  for  m  >  i,  we  can  write  any  polynomial 
of  degree  m  -- 1  as  a  linetir  combination  of  {p,(jc)  jlH'o  : 

M) 

Since  -  iajkj)p„(x) «  q„.i(x)  is  a  polynomial  of  degree  m  -  1  the  induction  hy¬ 
pothesis  implies  there  is  a  representation 


(.  m-t 

)  Pmix)  -  X 

Now  set  —  (ajkj  and  the  result  (3)  follows.  (Ref.  4:  p.  33] 

Using  the  theory  of  Fourier  series  developed  earlier,  we  next  determine  the  coeffi¬ 
cients  a,^  explicitly.  For  /  =  0, 1, ... ,  m,  let  c,^  «  and  let 

Then  by  construction,  { is  an  orthonormal  set  of  polynomials.  Writing  (3)  as 

m 

/■*0 

we  see  that  the  right-hand  side  can  be  interpreted  as  a  (terminating)  Fourier  expansion 
of  ^„(jc).  Hence  the  results  of  Example  15  in  Chapter  11,  Section  B  may  be  applied.  In 
particular,  by  (7)  in  that  section,  the  coefficients  arc  given  by 

”  <  9m .  *  f  9mW  0/W  »<Jf)  dx. 


Changing  back  to  the  old  variables, 
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-  */  <  <7m  .  J  ^mW  Mx)  dx. 


We  are  now  in  a  position  to  show  that  the  orthogonality  property 


(4) 


I 


can  be  expressed  equivalently  as, 


(5) 


b 

p^j(x)x'"w{x)dx^O ,  m<n. 


To  see  (5),  substitute  the  form  (1)  for  p,(jc)  into  (4)  where  m</i.  The  linearity  of  the 
integral  gives  (5).  On  the  other  hand,  since  x*  is  a  polynomial,  we  can  write  x**  as  a 
linear  combination  of  tlie  orthogonal  polynomials,  so  (5)  gives  (4).  Note  that  (5)  implies 
each  p„(x)  is  orthogonal  to  every  polynomial  of  lower  degree.  [Ref.  4:  pp.  33*34] 

B.  THREE-TERM  RECURRENCE  RELATION 

The  three^term  recurrence  relation  is  a  useful  result  which  holds  for  any  three  con¬ 
secutive  orthogonal  polynomials: 

(6)  p„(x) » {A„  X  +  R„K_i(af)  -  C„  p„^i{x)  , «  -  2, 3, 4, ... 


where  and  Q  are  constants  given  by 


(  \  ^  ^n-2  ^n-2  ^ 


The  recurrence  relation  is  valid  for  n  «  1  if  p_, «  0  with  Q  arbitrary.  In  this  case,  the 
formula  for  A„  also  holds  for  /r->  1.  (In  the  contrapositive  form,  this  statement  is  a 
powerful  tool  for  showing  that  a  polynomial  set  is  not  orthogonal.)  [Ref.  5:  p.  234] 

To  prove  this,  we  begin  by  considering  pj[x)  —  (A,/A,.|)  x /»,_,(x),  a  polynomial  of  de¬ 
gree  no  greater  than  (n—  1).  We  expand  it  in  terms  of  the  orthogonal  polynomials 
{pix)}U  via  the  technique  of  the  previous  section  to  obtain 
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(7) 


Pnix)  -  -T-®-  Jf  P„-|(jf)  -  2 


The  coeflicients  a,^i  are  determined  by 


<xp^\ix),pAx)>  ,  I  ^  a  -  1 


Because 


-J  xp^^ix)p^x)w{x)<lx^  <p^i(x),xp^x)>, 

it  follows  that  for  i  ^  —  3,  x  pj[x)  is  of  degree  no  greater  than  («  —  2).  Hence 

<Pn-i(x),xp^x)>  -  <xp^^(x)  ,p^x)>  •O  ,  I  i  a  -  3, 

since  orthogonal  to  every  polynomial  of  lesser  degree.  Thus  the  constants 

•••  f  2«ro,  leaving  o^i^,  and  With  this  knowledge,  (7) 

becomes 


Pnix)  -’r^xp^_y(x)  -  +  »^i^iP^xix)> 

'^n~l 

Setting  -  kjk^x ,  B,  -  ,  and  C,  -  -  a^,^,  then  rearranging  terms  gives  (6). 

To  determine  C,  explicitly,  we  write 


from  (8).  Since 


Q  -  -  <Pn-\ix)  ,Pn-.7(x)> 


XPn-2(x)  -  +  - 

^«<-2  r.  «-i  ,  1 

-T - +-J 


we  can  write 
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A 

Cii " Aw— 2  ^  ^1^11— iW  1(*)  A/jt-a pA'^)^  ' 


11-2 

z 

/-o 


Using  the  properties  of  the  inner  product  and  the  orthogonality  property  of  the 
polynomials,  we  obtain 


Cn 


An~a 

An-i  ^*-1 


Aw~2 

(*«)‘  V, 


as  given  in  (6).  [Ref.  6:  p.  8] 

Into  equation  (6)  we  substitute  the  expanded  forms  of  the  polynomials 
pMt  Pif-A*)  I  P»-A*)  (1)  and  the  constants  A,,  C,  from  above.  Equating  coelli- 

dents  of  gives  R..  Since  >  0  and  /t,  >  0,  it  follows  that  i4,  >  0  and  C,  >  0  and  the 
proof  is  complete. 

The  nonnegativity  of  the  constants  A,  and  C  is  important  for  the  converse  of  the 
result  in  (6).  Favard  showed  that  the  existence  of  a  three-term  recurrence  relation  in  the 
form  of  (6)  implies  that  the  polynomials  of  the  set  are  orthogonal  with  respect  to  some 
weight  function  over  some  interval  using  Stieltjes  integration  (Ref.  7]. 

We  observe  that  in  order  to  generate  the  polynomials  of  an  orthogonal  class  witli 
the  recurrence  relation  (6),  we  need  the  sequences  of  constants  and  C,  together 

with  two  of  three  consecutive  polynomials  in  the  class.  Other  techniques  provide  what 
some  authors  call  a  pure  recurrence  relation  requiring  only  two  of  three  consecutive 
polynomials  to  define  the  class,  because  the  constants  as  functions  of  n  are  contained 
explidtly  in  the  recurrence  relation.  The  recurrence  relation  derived  in  Chapter  I,  Sec¬ 
tion  A<1  for  the  Chebyshev  polynomials  is  an  example. 

C.  CHRISTOFFEL-DARBOUX  FORMULA' 

The  ChristolTel'Darboux  formula  is  an  important  identity  v'hich  can  be  derived  from 
the  three-term  recurrence  relation.  The  identity  is 


To  prove  (9),  note  that  from  (6)  we  have 

Py+, (x)  -  (/ly+,x  +  Bj+i)p/^x)  -  Cy+,Py_,(x) 
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which  we  rearrange  to  give 


(10) 

and  similarly 

These  recurrence  relations  are  valid  for  j  ->  0  if  we  set  C,  -  0  and  p.^(x)  -  p.^iy)  -  0. 

Multiply  (10)  through  by  p/^y)  and  (11)  by /i/x),  subtract  the  results,  and  then  mul¬ 
tiply  through  by  h,  to  obtain 

f»j  (jf -y)  pM)  p/y)  -  ^  [ Pj^iix)p/,y)  -  /^y+,(y)/»/x)] 

h  k 

+ [  Pj^ii*)p/.y)  -  /»/-i(^)/>/jf)]. 

Summing  overy  from  0  to  «  yields  a  telescoping  series 

'-y)^^PAx)Pj(y)  -  *7^  [  P^iix)Pniy)  -  P«(jf)/»^,(>')] 

JmO 

from  which  the  identity  (9)  follows.  (Ref.  4:  p.  39) 

Now  subtract  and  add  the  quantity  p,(jr)p^,(af)  to  the  numerator  of  the  right-hand 
side  of  (9)  and  let;;  tend  to  jc  to  obtain  a  limiting  case  of  the  Christofrel-Darboux  for¬ 
mula: 


2 ^[pM^ r " i^)Pnix) -y»\(jf)/>«+|(x) ]. 

>0  '^"+1 

We  will  use  (12)  in  the  next  section. 

D.  ZEROS  OF  ORTHOGONAL  POLYNOMIALS 

In  Chapter  I,  we  observed  that  the  zeros  of  the  Chebyshev  polynomials  {ri-r)}  are 
real,  distinct,  and  lie  in  the  interval  (—1,1).  This  principle  extends  to  any  class  of 
polynomials  orthogonal  on  tlie  real  interval  [a,6]. 
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To  SCO  this,  choose  n  >  0  tnd  suppose  that  p^x)  is  of  constant  sign  in  Then 
<  ,  ^jc)  >  #0,  which  contradicts  the  assumed  orthogonality.  Thus  by  continuity 

(and  the  Intermediate  Value  Theorem)  there  exists  a  zero  x,  •  ia,b). 

Suppose  that  x,  is  a  double  root.  Then  p,(jf)/(jf  -  Jfi)*  would  be  a  polynomial  of  de¬ 
gree  (a  —  2)  and  so 

0-  <p»(x),p«(jf)/(x-x,)^>  -  <  I ,  (p„(jr)/(jf-Jfi))^>  >0 

which  is  a  contradiction.  Thus  the  zeros  are  simple. 

Now  suppose  that  pJlx)  has  exactly  y  zoos  X|,  Xt» ... » Xj  •  («,h).  Then 

p„(xKx  —  JfiK-*  “  *2)  •••  “  *2^  ••• 

where  q^j  (jc)  does  not  change  sign  in  (a,b)  and 

< p^(jr) ,  (jr  -  JTjXjf  —  ...(.»->))>  "  *•  ”*  “  -*2)* ...  (or  -  .Ty)  > . 

Since 


and 

< P«(Jf) .  (Jf  -  -  Jf2)  -  (Jf  -  ®  <  a, 

then  it  must  be  that  j^n.  The  Fundamental  Theorem  of  Algebra  precludesy  >  a  and  so 
we  conclude^  —  a.  (Ref.  5:  p.  236) 

Thus  all  the  zeros  of  p,(;c)  are  real,  simple,  and  lie  in  the  interval  (0,6),  and  so  may 
be  ordered 

«  <  <  ^2/1  <  <  ”■  <  <  *• 

The  interlacing  of  zeros  of  p,(jf)  and  p,*,(jc)  follows  from  the  ChristofTel-Darboux  for- 
mula.  Recalling  that  the  leading  cocmcient  A,  is  positive  for  all  p,(x),  then  (12)  gives 

(13)  p;+i  -00  < a: <00. 

Let  u  and  v  be  adjacent  zeros  of  p,(af).  Then 

(14)  /»'»P'n(v)<0 
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since  the  zeros  ere  simple.  At  these  zeros,  inequality  (13)  reduces  to 


and 


Multii^y  these  two  inequalities  together  with  (14)  to  conclude 

so  p^iix)  has  a  zero  between  each  pair  of  consecutive  zeros  orp,(jr). 

Now  let  denote  the  largest  zero  of  p,(jr).  Observing  that  pXx)  oo  as  jr  oo, 
we  must  have  p',(z;J  >  0,  and  so  by  ( 13) 

But  Pm.|(x)  oo  as  X  oo,  so  Ah.i(x)  must  have  a  zero  to  the  right  of  x^.  Similarly, 
p^iix)  must  have  a  zero  to  the  left  of  X|^,  the  smallest  zero  of  /t,(x).  Thus  all  /t+l  zeros 
ofpUk(x)  are  accounted  for  and  interlace  those  ofp,(x)  [Ref.  8): 

*  ^  ^  ^  ^  ***  ^  ^  ^ ^  ^  ^  ^ 

£.  GENERATING  FUNCTIONS 

The  function  F{x,t)  having  a  formal  power  Mries  expansion  in  t 


ItmO 

is  said  to  be  a  ^entrating  function  for  the  set  (/Kx)}  (Ref.  9:  p.  129].  For  an  appropriately 
chosen  generating  function,  the  generated  set  {^(x)}  is  a  class  of  orthogonal  polynomials. 
By  defining  a  class  in  this  way,  properties  of  the  polynomials  can  be  derived  from  the 
generating  function  itself. 

For  example,  consider 


(15) 


Hxj)  -  (1  -  2jc»  +  -  £/.«(■ 

nvO 


28 


which  for  fixed  x  is  the  Taylor  series  of  f{x^)  centered  at  r «  0.  If  we  restrict  x  to  the 
interval  [  -1»1],  then  by  considering  the  singular  points  of  we  ntay  conclude  that 
the  series  is  convergent  for  |r|  <  1  {Ref.  4:  p.  28].  Thus  we  can  determine  X,{x)  for 
/i-0, 1,2,...  by 

(i«) 

Equation  (16)  yields  1,  /Kx)  — etc.  We  also  note  from 

(16)  that 

-1.’ 


We  now  derive  some  basic  properties  of  the/|(x)  from  the  generating  function  in  (15). 
1,  Recurrence  Relation 

Oiffercntiating  (15),  we  find 


(17) 

nm\ 

oe 

(18)  ^~(x-i)(\-2xi  +  iY'^-J^«Ux)i"-'. 

n^\ 

Since  (x  - 1)  4^  —  1 4^  "  have 

ox  or 

oo  oo 

tx  -  o^/.' W'"  -  ‘Y,  '•/.W'"'  - » 

n»l 

which  becomes 


f;  x/.-(x)«'  -  2  »Mx)t’  -  Ya'mi’*'  ■ 


Hal 


Huai 
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Since  ^x) «  1,  we  can  start  the  sum  on  the  right  side  at  x  *  C,  then  re-index  so  that  it 
starts  at  R  a  1  again.  We  then  find 

^  ** 

R^t  Rat 

Equating  coefficients  of  r*  then  gives 

(19)  (jf) .  »»^i. 

Rewrite  (17)  as 

OO 

(20)  (l-2xr  +  rV'^-2//(x)f"“‘,  t^O. 

Substituting  the  appropriate  expressions  from  (20),  (18),  and  (15),  respectively,  into  the 
identity 

(1  -  /’)(!  -  2xr  +  r^)“’^  -  (2r)  (x  -  r)(l  -  2xr  -h  r  V'*  -  (l  -  2xf  + 
now  gives 

which,  upon  rearranging,  bwomes 

CK»  oe  OP  OO 

W'"  -  Z 

Ri*l  R«l  Rail  RbI 

Equating  coefficients  of  r"  and  gathering  terms  yields 

(21)  (2r  +  I)/,(x)  -/,V,  (X)  (X) ,  R  ^  1. 

Substituting  (19)  into  (21)  gives 

x/«V)-/,V,(x)-(R+l)/,(x) 
which  by  a  shift  of  index  from  r  -»  r  —  1  becomes 


OO  OO 

-(2»X 

R«1  RaO 
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(22) 


("*)  1(^)  *  w  i  2. 


Substituting  again  from  (19)  gives 

(23)  (jf*  -  l)/«'(x)  -  nx/^ix)  -  »/«_,(x). 

Multiplying  (19)  through  by  (jt*  -  1),  we  have 

-  0/i'(*)  “  “  0^(*)  •  (**  ~  (-"f)* 

Now  substitute  for  (;f*  -  l)/l'(x)  and  (x*  ~  1)/;^,  (x)  from  (23)  and  gather  terms  to  get 

(24)  nMx)  -  (2rt  -  1)  x/^,(x)  -  («  -  l)/^a(x) ,  «  i  2. 

Equation  (24)  is  a  thrte-ttrm  recurrence  relation  for  {/|(x)}.  The  advantage  of 
this  form  is  that  beginning  with  f^x)  and  fix),  we  can  now  generate  any  member  of 
{/l('’f)}  by  iterating  (24)  and  thus  avoid  the  differentiation  in  (16).  Note  that  since 
A{x) "  1  •  yi(jf)  ■  (24)  implies  that  {iS(x)}  is  a  set  of  polynomials.  (Ref.  9:  pp.  159-1601 

2.  Ordinary  Differential  Equation 

We  continue  the  same  line  of  reasoning  to  extract  additional  irtformation  about 
(/;(x)}.  Differentiating  (22)  yields 

(25)  xy;.,  (X)  -//(X)  -  («  +  !)/;_,  (X). 

From  (19)  we  have  /Hi  (x)  and,  after  differentiating,  /Hi  (x).  Substituting  these  ex¬ 
pressions  into  (25)  gives 

-  («  -  l)/,'(x)]  -/"(x)  -  («  +  l)Cx/'(x)  -  /»/(x)] 
which  when  rearranged  becomes  a  second  order  ordinary  differential  equation 

(26)  (1  -  x')/,'(jc)  -  2jr//(jc)  +  ii(n  +  DAM  -  0. 

The  {/(x)}  are  solutions  of  (26)  for  « •»  0, 1, 2, ... .  [Ref.  9:  pp.  160- 161] 

3.  Orthogonality 
Rewiting  (26)  as 

(27)  ■*[(*-  +  "<»  +  ‘>/.W  -  » 
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we  recognize  the  structure  of  a  Sturm'Liouville  eigenvalue  problem.  Since  tlie  points 
*  ■  ±  1  are  singular  points,  we  require  that/^x)  andyj'W  be  finite  as ±  1.  From  the 
associated  theory  of  the  singular  Sturm-Liouville  problem,  we  conclude  that  the 
are  orthogonal  on  the  interval  [  -1,1]  with  weight  ftmction  M<jf)  -  1. 

Alternately,  we  combine  (26)  with 

■^[(1  +"('"+  i)/.W-o 


where  n^m,  to  obtain 

(281  [('  -''*)/«•«] 

+  W"  + 1)  “ 

Since 

[(I  -  *’){/-W/.'W  -/-MAWl] 

-/,W  x  [(1  -  -/.w  [(1  - 

we  can  write  (28)  as 

•^  [(>  -  *’){/»W/.'W  -/«'W/.W)]  +  +  «  -  m]f.(x)/Jx)  -  0 


(«  -  m)(<i  +  m  +  l)/»/„(;c)  -  [(l  -  J'’)I/«'W/.W  -A.W/«'Wl]- 

Integrating  from  x  -  -1  to  x  »  I,  we  obtain 

(/I  —  fn)in  +  m  +  1)  f  ”  •*  ){>An 

J-I 

Since  (1  -  X*) »  0  at  X  -  ±  I,  we  have 
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in  -  m)in  +  m  +  I)  f  f„ix)f„ix)dx  *  0. 
•'-1 

Recalling  that  conclude 


f  fnix)fmix)dx^0 

•'-1 

i.e.,  the  polynomials  (/Kx)}  are  orthogonal  on  the  interval  [-1,1]  with  respect  to  the 
weight  function  w(x)  =  1.  [Ref.  9:  pp.  173-174] 

These  are  the  Legendre  polynomials  that  were  previously  defined  using  the 
Gram-Schmidt  process.  Thus  they;(jf)  of  (15)  are  in  fact  PJ(x)  and  (15)  may  be  written 


(29)  = 

n-O 


establishing  the  equivalence  of  the  generating  function  definition  with  the  Gram-Schmidt 
definition. 

Legendre  and  Laplace  concluded  that  the  P,(x)  in  (29)  were  polynomials  of  de¬ 
gree  n  in  the  variable  x  by  examining  a  series  expansion  of  the  function 


(-ir(i/a.-„ ,, 

m!  (n-2m)!  ^  ’ 


They  reasoned  the  orthogonality  directly.  From  (29)  we  write 


(l-2;cr-f-r^)~'/^  =  ^  P^x)  r‘ 
and 


(l— 2x5+5^)“'^^  =  ^  Pj{x)  y. 

y=0 

Multiplying  these  power  series  together  via  the  Cauchy  product  formula  yields 
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f  2  /  2  '  z  -I 

I— 2xr+r  v/1— Zm+5  n-o/t-o 


oo  n 


OO  p  • 

■I]  X!  ww  --v. 

n»0  L  /+y*n 


Integrating  from  -1  to  1  with  respect  to  x  gives 


p  _ dx 

•'->  J\-2xr  +  r^  ./T^ 


2x5  +  s 


■E  E  [r  wwrf^]''V. 

n«0  /+y«n  L 


Through  tedious  calculation,  the  leA-hand  side  of  this  expression  becomes 

1  ^  + 

) —  *®8  ,  j —  » 

^rs  1  —  .^rs 

a  function  of  the  product  rs.  From  this  we  conclude  that  the  coefficients  of  the  terms 
in  the  series  on  the  right-hand  side  are  zero  when  /  and 7  difier,  i.e., 


Pix)PJix) 

•’-1 


dx  =  0,  i^j 


and  the  orthogonality  is  established.  Finally, 


or  P,(l)  =  1  and  so  the  P^x)  are  in  fact  Legendre  polynomials.  [Ref.  8} 

In  a  similar  fashion,  the  norm  of  the  Legendre  polynomials  can  also  be  obtained 
from  the  generating  function.  Let 


=  V'=/' 


Square  both  sides  of  (29)  and  integrate  with  respect  to  .r  on  [  —1,1] 
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P,W  r  dx. 

}  J-1  1  -  2jcr  +  r 

By  orthogonality,  the  left-hand  side  is  zero  unless  rn  *»  n,  while  the  right-hand  side  is 
integrable  in  closed  form: 


nmO  '  ’ 


This  function  can  be  expressed  as  a  difference  of  two  logarithms,  each  of  which  has  a 
convergent  Maclaurin  series  expansion  in  ( —1,1).  When  combined  this  yields 


Z 


naO 


C„t 


2n 


2 

2n+  1 


Comparing  coefficients  of  /*"  on  both  sides  gives  c,  =  . 

Since  there  is  no  systematic  theory  for  determining  generating  functions,  finding 
one  for  a  polynomial  class  can  be  a  problem.  The  work  above  bears  this  out.  Unfor¬ 
tunately,  the  proofs  above  do  not  easily  generalize  to  related  classes.  With  this  in  mind, 
let  us  summarize  the  key  steps  in  proving  that  (15)  is  a  generating  function  for  the 
Legendre  polynomials.  First  we  established  that /;(!)=«  1  for  n^O.  Next,  we  showed 
that 


for  /„(x)  a  polynomial  of  degree  n  in  the  variable  x.  Finally,  we  showed  that  these 
polynomials  were  orthogonal  on  the  interval  [-1,1] 

f  /nW/mW  =  0 ,  m¥>n. 

•'-1 

These  three  points  are  sufficient  to  show  that  the  generated  class  is  the  Legendre  class 
of  orthogonal  polynomials. 
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We  now  provide  an  ahemate  proof  attributed  to  Hcrmite.  Our  interest  in  this 
proof  is  mainly  the  technique  which  suggests  a  method  of  generalization  that  we  will 
take  advantage  of  in  the  next  chapter. 

We  begin  with  (15)  reproduced  below 


naiO 


Multiply  both  sides  of  this  equation  by  x*and  integrate  from  •!  to  1  to  obtain 


(30) 


"ij  \  —  Ixt  +  t  n»0 


Now  change  variables  from  or  to  >>  via 

(l-2jrr  +  /^)-’^^-l-0^ 


giving 


f(l-/) 

- - -  4.^^  dx  ^  -  ty)dy. 


The  lefl'hand  side  of  (30)  becomes 


I 


I 

-I 


(1-(K) 


{\-ty)dy 


or  simply 


Expanding  the  integrand  of  this  last  expression  via  the  Binomial  Theorem,  we  obtain 


(1-/) 


.k-J 


dy. 
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When  written  as 


y-o 


we  can  identify  the  form  as  a  polynomial  of  degree  ^  in  the  variable  t.  Comparing  this 
form  to  the  right>hand  side  of  (30),  we  conclude  that 


f  “  0 1  for/»>A:, 

i.  is  orthogonal  to  all  polynomials  of  lower  degree.  The  same  argument  as  before 
gives  7^(1) «  1  for  n  ^  0  and  the  proof  is  complete. 

F.  HYPERGEOMETRIC  SERIES 

The  term  "hypergeometric*  was  used  in  1655  to  distinguish  a  series  that  was  "be¬ 
yond"  the  ordinary  geometric  series  1  4-  x  +  .  In  1812,  Gauss  presented  the  power 

scries 


1  ,  ah  X  1)^(^  +  1)  +  2) 

c  1!  c(c+l)  2!'^  c(c+l)(c  +  2) 

c  0,  -1,  —2, ...  which  is  known  as  Gauss'  scries  or  the  ordinary  hypergeometric  series 
[Ref.  10]. 

Convergence  of  this  series  for  |x|  <  1  follows  directly  from  the  Ratio  test.  By 
Raabe's  test,  convergence  can  be  shown  for  |x|  »  1  when  (c  —  a  -  6)  >  0  [Ref.  1 1:  p.  5). 
Gauss  also  introduced  the  notation  2pii<^fh  \  c ;  x2  for  this  series.  Note  that 
]F,[a,6 ;  c  :x]  may  be  considered  as  much  a  function  of  four  variables  as  a  series  in  x. 
[Ref.  12:  p.  1] 

With  the  shifted  factorial,  the  ordinary  hypergeometric  series  can  be  expressed 


2fi[a,6;c;x] 


oo 
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Below  are  some  examples  of  important  Amctions  which  can  be  expressed  as  ordinary 
hypergeometric  series. 

Example  1:  log(l  +x)-jf ,/■,[!,  1  ;2;  -x] 

Example  2:  sin"*(x)  -  x jF,[l/2, 1/2 ;  3/2 ;  x*] 

Example  3:  tan-'(x)  -  x  jf,Cl/2, 1 ;  3/2 ;  -x*] 

The  generedized  hypergeometrie  series  is  formed  by  extending  the  number  of  pa¬ 
rameters,  an  idea  attributed  to  Clausen  [Ref.  1 1:  p.  40). 


OO 

,F,Ca„  aj, ... ,  a, ;  *2.  ••• . 


(^t)n(^)<i  •••  (4r)« 
•••  ' 


Note  that  since  (a)^,/(a),  >  a  +  a,  a  hypergeometric  scries  £  c;jr  is  characterized  by  the 
fact  that  the  ratio  c^Jc^  of  coefllcients  is  a  quotient  of  two  polynomials  in  the  index  n, 
i.e.,  a  rational  function  of  n. 

The  Ratio  test  can  be  used  to  show  convergence  for  all  values  of x  when  r^s  and 
for  |jr|  <  1  when  r«s  +  1.  When  r>s+  1,  the  series  diverges  for  all  x#0  and  the 
function  is  defined  only  if  the  series  terminates.  The  series  terminates  when  one  or  more 
of  the  numerator  parameters  a^  is  zero  or  a  negative  integer  [Ref.  11:  p.  45].  This  is  an 
important  characteristic  of  the  hypergeometric  series  that  will  be  used  later.  A  power 
series  that  terminates  gives  a  polynomial  which  is  defined  for  all  x.  In  this  case,  the  pa¬ 
rameters  5| . b,  may  be  negative  integers  as  long  as  the  series  terminates  before  a  zero 

is  introduced  into  a  denominator  term. 

Examples  of  the  generalized  hypergeometric  series  include  familiar  functions  such 
as: 

Example  4:  (1  +  x)*  -  ,fg[  -a ;  —  ;  — x] 

Example  5:  e* » g^oC  -  ;  -  ;  x] 

Example  6:  sin  x  —  x  ,f,C  -  ;  3/2 ;  -  x*/43 
Example  7:  cos  x  » ,F,[  -  ;  1/2 ;  -  x*/4] 

Example  8:  The  Bessel  function  of  order  a 

W2)*„f,[- :«+l: -*’/4] 

■'•W - riTTi) - 


where  dashes  indicate  the  absence  of  parameters,  i.e.,  when  r  =  0  or  s  =  0.  We  will  also 
use  a  common  alternate  notation 
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for  either  the  ordinary  or  generalized  hypergeometric  series.  [Ref.  12:  p.  4] 
1.  Chu-Vandermonde  Sum 
The  Chu* Vandermonde  sum 


jf,C a ;c ;  1] - 

is  one  of  many  useful  summation  formulas.  Since  this  one  will  be  used  in  a  later  chapter, 
the  proof  is  provided  below. 

Basically,  this  is  a  consequence  of  the  General  Binomial  theorem 


OO 


Starting  with  the  identity  (1  *-  jcHl  -  x)-*  -  (1  -  x)-*-*,  expand  both  sides.  Us¬ 
ing  the  Cauchy  product  on  the  leR  side  and  the  General  Binomial  theorem  on  the  right, 
we  have 


iDn  ’ 


A 

where  q,  -  2 


(aUbU 
(i)*(iu  ■ 


Equating  coefficients  of  oc*,  we  have 


h  (l)*(I)n-*  "  (1). 


In  order  to  express  the  left  side  as  an  ordinary  hypergeomctric  series,  jF,,  mul¬ 
tiply  both  sides  by  (1),  and  use  the  identity  (lyfl),-*  -  ( •-1)*(  -«)*  to  obtain 


(31) 


E 

i^-O 


( -n)kia)k 

(D* 


(a  +  b)„. 


Next, 
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-  ( -lA  -irA*x* + 1)  ...(»+»-*- 1) 

-  ( -lA  -»x  - 1)  ...(-»-«+* + 1)[  11*  ] 


(-*-«+ Do¬ 
using  the  above  result,  equation  (31)  becomes 


(-«)*(«)* 


1 


*-0 


(1)*  (-6-/I  +  1)* 


-(-ly 


'  (-6-n+l),  • 


Le»  Cl 


—  /I  +  1,  and  substitute  to  get 
l]-(-l)' 


ic)n 


( ~c  +  a  —  H  +  1) ...  ( — c  +  a) 
ic)n 

(c  —  a  + 1»  -*  1) ...  (c  -  a) 


ic)n 


(C  - 
ic)n  ‘ 
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IV.  JACOBI  POLYNOMIALS  AND  SPECIAL  CASES 


This  chapter  focuses  on  the  orthogonal  class  known  as  Jacobi  polynomials, 

Until  recently,  the  classes  of  orthogonal  polynomials  considered  "classicar  were  usually 
given  to  be  Jacobi,  Gegenbauer  (also  called  ultraspherical),  Chebyshev  (of  first  and  sec* 
ond  kind),  Legendre  (also  called  spherical),  Laguerre,  and  Hermite.  The  Jacobi 
polynomials  hold  a  key  position  in  this  list  since  the  remaining  classes  can  be  viewed  as 
special  or  limiting  cases  of  this  class.  Today,  the  classical  orthogonal  polynomials  are 
taken  to  be  special  or  limiting  cases  of  either  of  two  very  general  orthogonal  classes 
known  as  the  Askey^Wilson  polynomials  and  the  q-Racah  polynomials,  between  which 
we  will  establish  a  formal  equivalence.  Because  of  their  complexity,  these  classes  are 
described  in  Chapter  VI  after  the  necessary  additional  theory  has  been  developed. 

The  results  derived  in  the  text  that  follows  are  arranged  in  tabular  form  by  class  at 
the  end  of  this  chapter. 

A.  JACOBI  POLYNOMIALS 
1.  Definition  /  Orthogonality 

The  Jj’cobi  polynomials  T^^{x)  are  generated  by  applying  the  orthogonalization 
step  of  the  Gram*Schmidt  process  to  the  standard  basis  {1,  jr,jr', ... }  of  -1,1],  with 
respect  to  the  weight  function  given  by  a  continuous  beta  distribution  on  [  —1,1] 

Mjc;*,  /?)-(! -af)*(l+x)^ 


for  «>-l,^>— l,i.e.. 


f '  lt%)  -xfil  ^xfdx-lht  «]-' 

•'-t 


The  Jacobi  polynomials  can  also  be  represented  by  hypergeometric  series 


(1) 


(« + r-«» 


n  +  a  +  P  +  I  I— X 
flt+l 


This  set  of  polynomials  is  thus  standardized  (as  was  done  for  the  Legendre  class  in 
Chapter  II,  Section  B) : 
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We  shall  demonstrate  that  these  two  characterizations  of  the  Jacobi  polynomials  are 
indeed  equivalent. 

To  show  that  the  polynomials  defined  in  (1)  are  in  fact  orthogonal  with  respect 
to  M>(jr ;  a,  /8)  -  (1  -  ;f)*(l  +  xf  on  [  -1,1],  it  suffices  to  show  that  l^\x)  is  orthogonal 
to  one  polynomial  of  degree  m  for  m  »  0, 1, 2, 1.  This  is  because  any 
0  ^  m  ^  n  —  1,  can  be  expressed  as  a  linear  combination  of  such  polynomials.  While  any 
polynomial  of  degree  m  could  be  used,  we  choose  (1  +  jr)*  for  reasons  that  will  become 
apparent. 

To  establish  orthogonality,  we  consider 

<  .  (1  +  xf  >  -  f  ‘  (1  +  xf  (1  -  xf  (1  +  xf  dx. 

*'-1 

By  (1),  the  right*hand  side  becomes 


Lul  («+l)*A!2*  J-i  J 

kmO 

Using  the  last  result  of  Chapter  II,  Section  C.l,  the  change  of  variable  x"  1  -  2t  yields 


+  y  f/  (->■),(>■+«  +  />+ 1)>  r(*+<.  +  i)r(w+/t  +  i)  1 

I^LV  *!  /  («  +  l)»r(*  +  m  +  «  +  /!  +  2)  J 


Identities  for  the  gamma  function  allow  us  to  simplify  this  expression  to  the  form 
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r(m  +  ot  +  /?  +  2)  ^  (m  +  «  +  /?  +  2)*  At! 

*-0 


which  can  then  be  written 

+  I)  r~n,n-¥ci  +  fi+l  I 
r{m  +  a  +  P  +  2)  m  +  a  +  jJ  +  2  ’  J* 

The  Chu-Vandermonde  sum  Trom  Chapter  II  allows  us  to  write  this  as 

2*^^**^'r(a  +  >i+l)r(m +  /?+!)  (m+l-n)^ 

r(/fi  +  ot  +  /?  +  2)  (/ft  +  ot  +  /?  +  2)„ 


Thus, 


<  /»<*'^(JC),(1+Jf)"  > 


c„ ,  m  »  « 

0,  m-0, 1, 2, 


n-1 


where 


n\  2r***P^^  r’(n  -t-  g  -H)  r(n  +  /?  -f  1) 

"  r(2«  +  a  +  /I  +  2) 

which  justifies  the  orthogonality.  It  is  possible  to  extract  the  value  of  by  modi¬ 

fying  this  argument,  but  we  defer  this  computation  till  the  next  section,  where  it  will  be 
easier. 

2.  Ordinary  Differential  Equation  /  Rodrigues'  Formula  /  Norm 

We  begin  deriving  the  ordinary  differential  equation  for  the  Jacobi  polynomials 
by  noting  the  general  formal  result: 


(2) 


d  Oj,  ...  ,  O),  O)  ...  Qf  flj  +  1 ,  ...  ,  Q,.  +  1 

*  '''[i . 6,'^J"  ''■U  +  1 . i.+  l'"'.' 


This  can  be  seen  from  the  definition  of  the  generalized  hypergcoinetric  series  by 

writing 
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. ^  r1^  **  V  Jf* 

. 6/'J  dx  ^  {b,),...{b,),  k\ 

y  (<*\)k  >»  iOr)k  X*~* 

E**  ^**t^*-H  “•  JC* 

(^iWi  •••  (^xWi  ^• 


Noting  that  (a)*^,  -  a  (a  +  1)»,  we  have 


J_  prai,...,ar  1  ai...<v  y  (a,  +  Qj^ ...  (a^  +  l)it  ;r* 

<*c  '^'  k . ».  ■  'J  ■  »i  •■• ».  a  +  '>*  ■•■  <*' + *>*  *' 

from  which  (2)  follows.  Differentiation  of  the  hypergeometric  scries  is  justified  by  re* 
calling  that  the  power  series  is  convergent  for  |  x|  <  1  when  r  *  s  +  1. 

We  apply  this  result  to  the  Jacobi  polynomials  to  obtain 


d  ( 

(«+i).  .r 

—a,  a  +  a  +  /?  +  1 

1  -X  1\ 

'dx  ( 

a-Hl  ’ 

2  Jy 

/  -1 

\  («+i). 

-a  (a  +  a  +  j?  +  1) 

f— a  +  l,a  +  a  +  /?  +  2 

"(  2 

a+1 

’^■L  .  +  2  = 

(a  -h  «  +  /?  +  1)  (g  +  2)^_|  (»-  1)1  A 

2(a-l)!  («  +  2U 


l-x 

2 


which  simplifies  to 


(3) 


With  these  results  established,  we  now  consider  the  orthogonality  property  of 
the  class:  set 


(4) 


P^-%)  /■?'"(-»)  (I  +  x/*-  [*?• <!„,. 


By  (3),  we  can  write 
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'S:!'- m+l+it  f  -*)■(! 

•'-i 

Integration  by  parts  }ields 

[,  ["£  ('-»)■(>+*)')] 

noting  that  the  boundary  terms  vanish,  so  that 

J <’  ■  *>-’  +  *>'■'  * 

where 

(7)  it%)  -  <i(l+x)lt%)  +  />(l-*)/>^«W 

is  a  polynomial  of  degree  («  +  1).  Thus  by  (6)  we  may  express  as  a  linear  combi¬ 
nation  of  Jacobi  polynomials 


«+i 

;-o 

We  would  now  like  to  show  that  c,  -  c,  -  —  -  c,  -  0,  so  that  only  the  last  term 
of  the  summation  survives.  The  constants  ,  y  —  0, 1,  2, ... ,  n  can  be  determined  by 
substituting  this  expression  for  into  (6)  and  using  (4). 

4';?-/'  if-%) .'$''(x)(i ^xfdx 
'*-1 

'•+1 

-„+^+p  Z'J ,  '’■“we  -xr'  ii+xf-' rfx. 

For  each  m  s  n  -  1,  by  (4)  the  left-hand  side  7*^  »  0,  but  the  right-hand  side  is 
zero  unless  j^m  +  1  ^  Thus 
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which  means  from  (6) 


^  r 

m+a+fi  "**  '•H-i.w-i 
and 
(») 

To  determine  the  remaining  constant  <^1,  we  note  that  by  (7)  and  (9) 

(i-x’)  - « (i+i)  it%) + f  (1-*)  f^-%)  -  c„,  />*:;’■  '-'W 

Letting  x«  1,  then  from  the  hypergeometric  series  deOnition  of  the  Jacobi  polynomials 
we  have  ^(1)  ■  which  gives 

It! 

-2a  -c^, 

and  so 

(10)  c„+,--2(«+l). 

Combining  (5),  (6),  (9),  and  (10),  we  have 

(11)  [/><••«(:,) (!-*)■  (l+x)'’]- -2 (n+l)[/>JJ;'’'’-'’(jt)(I-x)*-'(lW']. 

To  obtain  the  second  order  ordinary  difTerential  equation,  change  rt-*/i-l,  a-*a+l, 
and  P  -*  P+l  in  (1 1)  to  give 

(I-Jf)*"^'  (l+;c)'^'^'  ]  -  -2  «  (l-;c)“ (I+.r)''. 

Then  by  (3), 
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[  „+.+/)+ 1  ]  "  -2  "  ('-•')“  d+J')'. 

i.  e..^*  P^n'^\x)  satisfies 

■£  r  (!-:<)■■"  (l+x/+'  ]--n(«+K+/!+l)(l-Ar)'‘(l+;t)'y. 

Using  the  product  rule  to  expand  the  left-hand  side,  we  have 

(1-^:)“  (l+Jc)'"  [-(a+1)  (l+:c)y  +  (A+D  (l-Jf)/  +  J 

=  -n  {n+a-^p+l)  (l+x/j^ 

which,  when  simplified,  becomes  the  second  order  ordinary  differential  equation  for  the 
Jacobi  polynomials  j»  =  : 

(12)  ( 1  -Jf j'"  +  [  {0-0.)  -  (a+0+2)  jc  ]  y  +  « {n+a+P+ 1 )  >»  =  0, 

The  reader  is  invited  to  compare  this  result  with  the  second  order  ordinary  differential 
equation  for  the  Chebyshev  polynomials  {r,(jc)}  given  by  Equation  (5)  in  Chapter  I. 

By  iterating  ( 1 1)  A:  times,  we  obtain 

-4-  [  /’?■'”«  (1-Jr)“  (l+Jf)^  ]  -  (-1)*  2‘  (<1+1),  [  Ptt’  (1-^)-*  (l+4^-‘  ]. 

ax 

Setting  a  -*  a+A,  and  0  -*  0+k  gives 

2  {n+l)k  dx'^ 


or  equivalently  in  terms  of  the  weight  function, 

w(.v ;  a,  0)  P^:ii\x)  =  ^  ^  I 

2  (w+I)t  dx 


the  general  Rodrigues  formula  for  the  Jacobi  polynomials.  Letting  /» =  0  and  noting 
K'  “Xx)  =  1  gives 


(l-;c)^(l+.^)''7’r 


,P  p('*.  P)( 


(-If  /  r,, 


{l+xf 
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or 


w(^ ;  a,  /?)  Pf’  -77-  [  vv(jr ;  a+/c,  P+k)  ], 

k\  dx^ 

the  classical  Rodrigues’  formula  for  the  Jacobi  polynomials. 

Finally,  we  can  use  these  ideas  to  obtain  the  value  of  By  direct  com¬ 

putation, 


•'-1 

_  r(a-n)  r(/?-n) 
r(a+^+2) 


We  will  use  this  below.  Combining  (8),  (10),  and  (4)  when  m  *  /i  yields 


n+a+p  L  '»+•  -I  • 


Making  the  changes  n  -+  «-l,  a  -*  a+1,  and  P  -*  jS+I,  this  may  be  rewritten  as 


Iterating  this  relation  n  times  produces 


^(2w-t-«-^/?+l)  r(/t+a+l)  rjn+p+i) 

2^"  «!  r{n+a+p+l)  r{2n+a+P+2) 

2«+/?+i  r(/t4-tt-n)r(«+/?+i) 

(2/i+a+/J+l) «!  r(/i-l-a+/?+l) 


3.  Generating  Function 

As  was  mentioned  in  Chapter  III,  finding  a  generating  function  for  a  class  of 
orthogonal  polynomials  can  be  a  challenging  task.  Fortunately,  a  generating  function 
for  the  Jacobi  polynomials  can  be  found  by  mimicking  a  technique  attributed  to  Mcrmite 
in  his  work  with  the  Legendre  polynomials.  (See  Chapter  HI,  Section  E.3.) 

We  seek  a  generating  function 
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(13) 


n<BO 

for  constants  c„,  /t  =  0, 1, 2, ,  with 

f  ‘  x'‘il-x)‘‘il+xfdx^0,k<  n. 

•'-I 

To  begin,  we  multiply  (13)  by  x*  (1  -  jr)*  (1  4-  and  integrate  from  -1  to  1  to  obtain 

fc 

(14)J’^  x'^FMd  -jc)“(l  +4r)'’</j  =  ^c,j^  +;c)'’<fa  j  t". 

Note  that  the  summation  of  the  right-hand  side  is  from  0  to  A  because  the  orthogonality 
property  makes  each  term  zero  for  «  >  A.  We  observe  that  the  unknowns  to  be  deter¬ 
mined  in  (14)  are  F(r,x)  and  c„.  When  one  is  given,  the  other  can  be  found,  so  we  first 
consider  the  left-hand  side  of  (14) 


(*'  x*/'(r,x)(l-x)“(l-fx)^<Af. 
Setting  ./ 1— Ivr+H  m  l—ry  and  substituting  for  jf 


f I  r  (1-/)  f  >  c-j'') <•  \  r  (•-/) ^  > r  c-/)  ^  > 

^('■•3^— 2— )[‘-3 - ['+■>'+--2— J 


(IV) 


(IV) 


ry)  dy. 


Factoring  out  the  terms  (1-^)*  and  (!+>'>',  we  have 


£,  (l-^r  (1+3)'’  f 

We  would  now  like  to  make  a  judicious  choice  of  F  in  order  to  facilitate  calcu¬ 
lations.  Following  the  lead  of  the  Legendre  polynomials  (Chapter  Ill,  Section  E.3) 
suppose  F  is  such  that 


49 


(16) 


then  the  integral  in  (15)  becomes 

We  note  that  the  integrand  is  a  polynomial  of  degree  k  in  the  variable  r  as  is 
the  right-hand  side  of  (14).  From  (16),  we  get 

which  becomes,  via  l-2Jfr-l-r*  =  (l-ry)^  the  generating  function  for  the  Jacobi 
polynomials; 

f  (r,;r)  =  2“+^  (l-2;cr-br^)'’‘/^  [l~r+(l-2;rr-l-r^)‘/^]"*  [Hr+(l-2xr+r^)‘'^]”^. 
Using  this  generating  function  in  (13)  gives 

2  c„  J'"  =  2“+^  {l-2xr+r^)-^'^  [l-r+(l-2.xr+r^)‘/^]"“  [l+r-l-(l-2.Tr-t-r^)'/2]"^ 

/»=0 

from  which  we  can  determine  c,  by  setting  x*  1  and  recalling  that  f<‘’^>(l)=  , 

nl 

Thus, 

2  [2~*  (1-r)"*]  2"^  -  (l-r)~“-‘ 

«-o 

which  by  the  Binomial  Theorem  becomes 

V  («+»«  «  V  («+>)«  „ 

nstO  /i=sO 

and  so  c,  =  1  for  all  n.  [Ref.  8] 
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B.  SPECIAL  AND  LIMITING  CASES 

With  the  structure  of  the  Jacobi  polynomials  established,  we  turn  now  to  the  role 
of  the  parameters  a  and  p. 

1.  Special  Cases 

For  certain  choices  of  the  parameters  a  and  P,  we  find  that  the  classes  previously 
examined  and  several  new  classes  are  produced  as  special  cases  of  the  Jacobi  class. 
These  subclasses  inherit  the  structure  of  the  parent  class  which  oflen  provides  a  direct 
way  to  establish  specific  properties  (i.e.,  polynomial  nature,  orthogonality,  etc.) 

Table  1  provides  the  choice  of  parameters  for  selected  classes. 


Table  1.  PARAMETERS  FOR  SELECTED  CLASSES 


Class 

Parameters 

Jacobi 

a>— 1,  P>—\ 

Gegenbauer 

a  =  /l  =  ;-l/2 

Chebyshev, 

First  Kind 

Chebyshev, 

Second  Kind 

II 

II 

« 

Legendre 

a^P  =  0 

2.  Limiting  Cases 

In  this  section,  we  briefly  examine  the  Laguerre  polynomials  and  the  Hermite 
polynomials.  Using  the  hypergeometric  series  definition  of  the  Jacobi  polynomials,  we 
show  how  the  Laguerre  class  is  a  limiting  case  of  the  Jacobi  class. 
a.  Laguerre  Polynontiais 

The  Laguerre  polynomials,  are  defined  in  terms  of  hypergcometric 

series 

(.7) 

The  first  relationship  to  establish  is 
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lim 


/J-»oo 


To  show  this,  we  begin  with  the  hypcrgeometric  series  representation  for  the  Jacobi 
polynomials 


2x 


When  y  »  1-  ,  we  have 


"  V  ArUj'‘L  «+i  ■ />  J 


Writing  out  the  power  series,  we  obtain 

2£.\  /«+«\Y  i-n)kin+<i+P+l)k  X* 

"  \  P  n  («+l)*A! 

_  /  n+a.  \  ^  {n+(i+fi+l)i( 

“V  «  /jS  («+l)fcA! 


In  the  limit  as  /?  -+  oo,  the  ratio 


{n+(x+p+l)fc 


and  so 


(a+I)*^! 


completing  the  proof.  [Ref  13:  p.  103] 

To  obtain  the  orthogonality  relation  for  the  Laguerre  polynomials,  we  start 
with  the  orthogonality  relation  for  the  Jacobi  polynomials 


•*-1 
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Letting 

Passing  the  constants  outside  of  the  integral  and  dividing  through  by  them  leaves  us  with 

J"  '”(1-^ )  e  «(i-  ^ )  X-  (i-  j-  y  A  -  0 . « # 

Taking  the  limit  as  P-*oo,  we  obtain  the  orthogonality  relation  for  the  Laguerre 
polynomials 


(18) 


(19) 


Jn 


b.  Hermite  Poiynomials 

The  Hermite  polynomials,  H^x),  are  defined 


H„ix)^{2x)\Fc 


„  -«/2,(-«+l)/2 


0 


m¥‘n. 


In  a  fashion  similar  to  that  for  the  Laguerre  polynomials,  the  Hermite  polynomials  are 
a  limiting  case  of  the  Jacobi  polynomials  via  the  Gegenbaucr  polynomials.  Specifically, 

H„ix)  =  «!  Urn 


which  allows  a  derivation  of  the  orthogonality  relation 


(20)  J  H„{x)  H„ix)  e~^  dx=^0 ,  m^^n 

”00 

from  that  of  the  Jacobi  class.  [Ref.  13:  p.  107) 

C.  DISCRETE  EXTENSIONS 

We  turn  now  to  orthogonal  polynomial  classes  which  use  the  discrete  inner  product 
introduced  in  Chapter  II,  Section  A.2 
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</.  g> 


AT 


XmO 


Here  instead  of  a  continuum  the  support  of  the  weight  function  is  concentrated  of  a  fi¬ 
nite  set  of  discrete  mass  points  {0, 1,...  .AT).  The  polynomial  classes  of  particular  in¬ 
terest  are  the  Hahn,  dual  Hahn,  and  the  Racah  polynomials. 

1.  Hahn  Polynomials 

The  Hahn  polynomials  -  actually  discovered  by  Chebyshev  -  were  independently 

% 

realized  by  physicists  working  in  angular  momentum  theory  via  3-j,  or  Clebsch-Gordon 
coefficients  [Ref.  8].  We  define  this  class  by  the  generalized  hypergeometric  series 


(21) 


Qnix’,a,  P,  N) 


— jr,  n-|-a+j5-f-l 
a+l,-Ar 


for  «  >  — 1 ,  fi>  —1,  where  N  is  &  positive  integer  and  n  =  0, 1, ... ,  N.  From  the  power 
series 


we  note  that  the  variable  jr  does  not  appear  where  we  have  come  to  expect.  Since 


(-jr)*-(— ;f)(— jc+l)  ...  (— x-f*/i-l) 
.(-!)*  W(*-l)  ...  (x-*+l), 


we  conclude  that  QJ[x ;  a,  AO  is  a  polynomial  of  degree  n  in  the  variable  x.  Because 
O^n^N,  this  set  of  ^+1  orthogonal  polynomials  is  finite  for  fixed  a  and  fi.  [Ref.  14: 
p.  331 

The  Hahn  polynomials  satisfy  the  discrete  orthogonality  relation 


AT 

(22)  ^  Qjx  p,  N)  Q„ix  ;a,  P,  N)wix',(x,  p,  N)^0 ,  m^n 

where  the  weight  function  is  given  by  a  hypergeometric  distribution 
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(23) 


(i)x  (l),v-x 


We  note  that  the  weight  function  may  also  be  written 


(24) 


/?,  N) 


ifi+lh  («+l)x(-iV)x 

(1);,  (I);c(“^-/»)x 


We  introduce  the  first  forward  difference  operator  acting  on  x 

-/n(^+l)  ~/,(jf) 

as  a  discrete  analogue  to  differentiation.  Since 


E 


A/,(/)  -/«(x+l)  -/«(0). 


(a  discrete  analogue  of  the  Fundamental  Theorem  of  Integral  Calculus)  the  first  forward 
difference  operator  is  a  discrete  inverse  of  the  summation  operator.  This  difference  op¬ 
erator  is  used  in  the  Rodrigues'  formula  for  the  Hahn  polynomials,  which  can  also  be 
written  in  terms  of  the  weight  function  wfjr ;  a,  N)  with  shifted  parameters  as  was 
done  for  the  Jacobi  polynomials. 

Two  limiting  cases  of  the  Hahn  polynomials  are  of  particular  interest.  In  the 
first  case,  replace  x  by  Nx  in  the  interval  of  orthogonality.  This  in  effect  places  the 
support  of  the  weight  function  on  the  ( iV+l)  equally-spaced  points  { 0, 1/zV,  2/fV, ... ,  1 } 
in  [0,1].  As  N-^oo,  the  set  of  discrete  mass  points  tends  towards  the  full  interval 
[0,1] ;  we  expect  that  this  structure  is  reflected  in  the  discrete  orthogonality  relation  (22) 
becoming  an  integral  orthogonality  with  respect  to  a  continuous  weight  function  on  that 
interval.  Rewriting  the  weight  function  in  (23)  as 


M.<jf  ;a,  p,  fV)  = 


r(.v+l+«) 

r(.r+l) 


r(  N-x+\-^p) 

n  iV-jc+i) 


and  using  a  consequence  of  Stirling's  Formula  [Ref.  15:  p.  257] 
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1, 


limr"* 


r(f+fl) 

r(/) 


we  see  that 

r(a+l)r()3+l)lim  0, 

A^-*oo 

which  we  recognize  as  the  continuous  beta  distribution  for  the  Jacobi  polynomials  nor¬ 
malized  on  [0,1]. 

This  is  indeed  the  case,  and  can  be  easily  verified  by  a  direct  computation  on  the 
hypergeometric  series  definitions  (21)  of  J  “•  A.  (0  of  given  in  Sec¬ 

tion  A.l,  i.e.. 


jima(iVx;«,  /!,  />?■ '"(l-2x). 

Thus  the  Hahn  polynomials  may  be  viewed  as  a  discrete  analogue  and  generalization  of 
the  Jacobi  polynomials.  [Ref.  14:  p.  36] 

The  second  liiniting  case  gives  rise  to  an  interesting  class  of  polynomials  which 
has  applications  in  coding  theory  [Ref.  16). 
a.  Krawti'houk  Polynomials 

For  0  <  p  <  1,  let  a  =  pr ,  —  (1— />)  t  in  (21),  then  take  the  limit  as  t  tends 

to  infinity  to  obtain  the  Krawtchouk  polynomials,  i.e., 

limC;,(ar \pt ,  (1-/?)  / ,  fV)  - 

t-*oo 

(25)  jf,  :  7  ]  -  A:.(x  ;  p.  AT ) 

[Ref.  14;  p.  38).  A  similar  limiting  argument  applied  to  the  weight  function  shows  that 
the  Krawtchouk  polynomials  are  orthogonal  with  respect  to  a  binomial  distribution: 

(26)  J]k„{x-,p,  N)K„ix;p,  /V)('  )/(l-p)'^'*-0 ,  m^^n. 

This  class  of  polynomials  possesses  an  inherent  synmietry,  the  structure  of 
which  can  be  generalized  to  form  other  orthogonal  polynomial  classes.  To  this  end,  we 
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turn  now  to  the  characteristic  of  duality.  For  suitably  defmed  functions  mW  and  v(j:),  two 
classes  of  orthogonal  polynomials  { pXtdx)) ),  and  { 9;(v(jr)) }  arc  said  to  be  dual  if 

P;c(A*(rt))  -  9n(v(jf)). 

That  is,  interchanging  the  roles  of  the  degree  n  and  the  discrete  variable  x  in  one  class 
produces  the  other.  The  Krawtchouk  polynomials  provide  an  example  of  an  orthogonal 
class  that  is  seff-dual,  i.e., 


N). 

This  is  clear  from  the  hypjrgeometric  series  definition  (25)  given  above.  It  therefore  may 
seem  reasonable  to  suspect  that  there  exists  a  class  of  orthogonal  polynomials  dual  to 
the  Hahn  polynomials.  Such  a  class  does  in  fact  exist,  and  it  is  this  dual  class  which  we 
next  examine.  (Ref.  17:  p.  657] 

2.  Dual  Hahn  Polynomials 

The  dual  Hahn  polynomials,  ^.(jc) ;  y,  5,  N ),  are  defined 


(27) 


/?„(  >l(jr) ;  y,  5,  N) 


-Of,  x+y+^+1 

y+l,-iV 


where  ,i(jr) ■=.i:(jf4-y+5+l),  and  «=>0, 1,...  These  polynomials  satisfy  the 
orthogonality  relation 


(28) 


s 

^  ; y, <5,  N)  Rn{X{x)\y,5,  N)M.x\y,5,  /V)  =  0,  m^n 


x-o 


where 


,  .  ,  „ .  (y+a+1),  i{y+d+3)l2):  (y+ 1),  j^N ).  (-1)" 

,2^  ’  ’  (l)x((y+«5+l)/2)^(.5+l),(y+J+.V+2)^ 

/.vx  (y+«^+i)x (y+  Ox (~^x  2jf4-y+5+i 

^  il)^(S+l)^{N+y+d+2)^  y+d+l  • 

Note  that  R„{  Aix) ;  y,  5,  iV)  is  a  polynomial  of  degree  n  in  the  "variable"  /1(j:).  The  rea¬ 
son  can  be  seen  directly  from  the  hypergeomctric  scries 


57 


58 


Itji(n);c^  f.  A  n 

This  duality  is  also  reflected  in  the  recurrence  relations  and  difference  equations  of  the 
two  classes.  By  interchanging  the  roles  of  n  with  x,  et  with  y,  and  p  with  6,  the  recur¬ 
rence  relation  for  one  class  leads  to  or  can  be  extracted  from  the  difference  equation  for 
the  other  class.  [Ref.  14;  p.  37] 

3.  Racah  Polynomials 

Both  the  Hahn  and  dual  Mahn  polynomials  can  be  unified  as  special  cases  of  a 
single  larger  class,  the  standard  notation  of  which  is  similar  to  that  of  the  dual  Hahn 
class.  The  Racah  polynomials,  RS.  /l(jr) ;  a,  /S,  y,  6\  are  defined  by  the  hypergeomctric 
series 


(30) 


P,Y,6)^aP: 


Fr-"’ 


-n, «+«+/?+ 1,  —X,  x-l-y-l-5+1 


a+1,  )?+i+l,y+l 


where  i(jr) » jc(x-l-y+<5+l),  »•  0, 1, ....  AT,  and  one  of  «+l,  )?+5+l,  or  y+1  equals  -N. 
(Physicists  understand  these  objects  via  6-j  symbols.) 

The  orthogonality  relation  for  the  Racah  polynomials  is 


ff 

(31)  2  R„(  ^(of) ; «,  p,  y,  a)  /?„(  i(jf) ;  a,  P,  y,d)w{x',a,  P,y,6)^0 ,  m^n 

x-o 


where 


(y+<S+  l)x  {iy+M)l2)x  (a+l)x  ( P+M)^  (y-f  1)^ 
"  (l),((y+<5+l)/2),(y+a-a+l)j,(y-/?+l),(a+l), 

(y+J+l);t  (tt+l)y(/?+a+l)x(y+l)x  2jc+y+a+l 

"  (i)jc(y+^-«+i)x(y-^+i);c(«5+i)x  y+s+i 


[Ref.  19:  p.  24).  Note  that  w(jr ;  a,  p,  y,  5)  has  a  "vrcll-poised"  structure.  This  means  that 
the  pairwise  sum  of  numerator  and  denominator  parameters  is  constant,  i.c., 
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(y+5+l)  +  (I)-y+a+2 

(a+1)  +  (y+5~«+l)  ■  y+i+2 
(/J+J+I)  +  (y-^+i)  -  y+5+2 
(y+1)  +  (^+1)  -  y+^+2. 


The  same  could  be  said  for  the  weight  Amction  of  the  dual  Hahn  polynomials.  (To  say 
that  a  hypergeometric  series  is  well-poistd  means  that  for 


F,h’ 4 

L  ••• » J 


the  parameters  a  and  b  satisfy  the  following  relation 

a|+l  ■«  Oi+bi  •  —  ■■ 


Well'poisedness  is  an  important  property  of  certain  summable  hypergeometric  series.) 

From  the  hypergeometric  series  definition  (30)  of  the  Racah  polynomials,  we 
find  that  when  the  roles  of  x  and  n.  a  and  y,  and  and  ^  are  all  interchanged,  the  series 
is  unchanged.  Thus,  like  the  Krawtchouk  polynomials,  the  Racah  polynomials  are  self¬ 
dual. 

To  recover  the  Hahn  polynomials  as  a  limiting  case  of  the  Racah  polynomials, 
let  y+1  ■  —N  and  6  ~*oo.  Thus  by  formulas  (30)  and  (21), 


lim  4^3 

0-*oe 


•rt,  n+a+^+l,  -jc,  jf+y+d+1 
#+l,  P+S+l,  — iV 


— rt+a+jJ+1,  — jf 
0+1,  -iV 


and  then  by  formulas  (32)  and  (24), 


lim  w(jf ;  0,  p,  y,  d) 
0^00 


(^+1)a' 


w(x Id,  P,  N ). 


Likewise,  by  letting  «+l  »■  —iW  and  /J  -♦  oo,  we  obtain  the  dual  Hahn  polynomials. 

Figures  5  and  6  together  with  Tables  2  and  3  provide  a  hierarchy  of  classes  dis¬ 
cussed  in  this  chapter.  Tables  4-15  summarize  information  about  selected  classes  in  this 
hierarchy  [Refs.  14,  12.  20,  21  ,  5,  13,  9 , 22j. 
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Figure  6. 


Jacobi 


Gegenbauer 


First  Kind  Second  Kind 


Hierarchy  of  Special  Cases 
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Table  2.  TABLE  OF  CLASSES 


Class 

Symbol 

Racah 

/?„( /l(x) ;  0,  P,Y,S) 

Dual  Hahn 

R„iXix)-,y,5,  N) 

Hahn 

Q„{x\(x,  p,  N) 

Jacobi 

pt-%) 

Krawtchouk 

K„ix‘,p,  N) 

Mcixner 

M„ix ;  /?,  c) 

Laguerre 

l'X) 

Charlicr 

C„{x ;  o) 

Hermite 

//nW 

I  lypergeometric  Series 


— rt,  — Jf,  X  +  0L  +  P+l.~ 

y+l,  -N 


r  r~'^*  w  +  a  +  /?  +  1 .  , ' 
0+1,  -/V 


(a  +  On 


Table  4.  RACAH  POLYNOMIALS _ 

Symbol:  /?,( ^(.r) ;  a,  /?,  y,  S) 

Interval:  For  N  a  positive  integer,  =  0, 1, ... ,  N. 

Weight: 

/  o  e\  (y+<S+l)jc  {{y+6+‘i)l2)x  (°^+0;c  ( (V+O;: 


Norm: 

ff 

^  [  R„i  ;.(jr) ;  a,  /?,  y,  <5)  f  M<jf ;  a,  /?,  y,  <5) 

x>0 

U  (^+°‘+/^+0n  (  A+On  («~^+i)n  (tt+/?~y+l)n 
(«+/?+ 2)2„  (a+l)„  (  ^+5+I)„  (y+l)rt 

where  if 
a+l=-iV, 


(y-l-<5-f2)/y(--/?)y 

iy-P+lhiM)ff  ' 

or  if 

/?+<5+l=-A^, 

Af_  (y+*^+^)vv(<^~“).v 

(y+5-a+l)^(<5+l);v  ’ 


or  if 

y+1  =-iV, 


Table  5,  RACAH  POLYNOMIALS  (CONTINUED) 


Hypergeometric  Series: 


R  ( Mr)  •  a  S  v  S)^  n+a+fi+l,  -x,  x+y+S+l , 

K„(  /AX) ,  a,  p,  Y,  ,  IJ 

where  i(jf)  =»  j:(jc+y+5+i),  «  «  0, 1, ... ,  iV, 
and  one  of  a+1,  )ff+5+l,  or  y+1  equals  -N. 

Recurrence  Relation: 

Xix)  R„i  Xix) ; «.  p,  y,  <5)  =  a(n)  [  /?„+,(  X{x) ; «,  /?,  y,  «5)-/?„(  ;.(;c) ;  a.  y.  S)  ] 

-c(n)  [  R„i  ^(x) ;  ot,  p,  y,  ;.(x) ;  a,  P,  y,  5)  ] 

where 

.A  _  (w+g+^+l)  (w+g+1)  (n+jg+^+1)  (w-fy+l) 

^  ^  (2/i+g+/?+l)  (2n+0L+P+2) 

,  V  n  (n+p)  (w+g+/?— y)  (/i4-g— 

(2/i+g+yJ+l)  (2w+a+y?) 

Difference  Equation: 

Xin)  R^i  2(n) ;  a,  p,  y,  S)  -  ^(x)  [  .l(«) ; «,  p,  y,  syR^^i  >!(«) ;  g,  p,  y,  6)  ] 

-C(x)  [  /?jj(  Xin) ;  g,  p,  y,  S)-^R^_^{  Xin) ;  g,  p,  y,  S)  ] 

where 

A  ^  (jf+y+^+1)  (jf+y+1)  (x+^-fjg+l)  (jr+g+1) 

^  ^  (2.v+y+«5+l)  (2x+y+5+2) 

.  X  (x+<$)  (x+y+<?-a)  jx+y-p) 

(2x+y+<5+l)  (2x+y+5) 


Table  6.  DUAL  HAHN  POLYNOMIALS 


Symbol:  /?,( AW  ;  y.  <5,  N) 

Interval:  For  ^  a  positive  integer, x  —  0,\,...,N. 


Weight: 


vi<x;y,5,  N) 


(v+s+ih  {iy+s+m)x  (y+i)x  i-^)x  (-1)' 

(l)jc  ((y+^+l)/2)jt  (^+i)jt  (y+5+iV+2);c 


Norm: 


2][«„(;,(Ar):y,i.  N)]^w(x-,,.S,  iV )  “  [  ("  J’’)  ( VtO  ]' 


Hypergeometric  Scries: 


where  A(x)  =  x  (jf+y+6+1) 

Recurrence  Relation: 

-AW  R„i  A W ;  y.  <5.  )  =  m  [  AW ;  y.  <5,  N  )-R„(  AW ;  y,  <5,  )  ] 

-D{n)  [  R„i  AW  ;y,s,  N  )-/?„_,(  AW ;  y.  <5,  iv )  ] 


where 


Bin)  *  ( N—n)  (y+l+zi) 
Din)  =  n(  N+l+d-n) 


Difference  Equation: 

-«/?;,( AW  ; y. <5,  N)^bix)[R„i^ix+l)]y,5,  N)-R„iAix)\y,d,  N)] 

-dix)[R„iAix);y,d,  fV)-R„(  A(.r-1) ; y, «5,  N)] 


where 


<f(x)  = 


(jc+y+^-t-l)  (jf+y+l)  (  N—x) 
(2jf+y+<54-l)  {2,T+y+5+2) 
jr(jf+^)(jc+y+0+iV+l) 
(2jf+y+«5)  (2jf+y+^+l) 


Table  7,  HAHN  POLYNOMIALS 


Symbol:  <2,{jr;a,  N) 

Interval:  For  N  a  positive  integer,  x  =  0, 1, ... ,  N. 


Weight: 


(A+Dv 

(l),v 


H)A-N-Ph 


Norm: 

OO 

^[Qni.x\Oi,  /?,  /?,  AT) 

x»0 

/iV+a+/?+l\  /Ar+a+/?+l+n 
a+/?+l  \  iV  A  n 

2/i+g+^+l 

r(g+i)r(g+;?-n)  r(/»+/?+i)  r(«+i) 

r(/?+i)  r(rt+a+i)r(rt+a+iS+i) 


Hypcrgeometric  Series:  For  g,  j8  >  -1  and 
n  =  0, 1, ... ,  N, 


au:.,  /I, 


Table  8.  HAHN  POLYNOMIALS  (CONTINUED) _ 

Recurrence  Relation: 

-xQ„(x;a,  fi,  /V )  =  A(/i) [ C„+i(jr ; a,  /?,  N)-Q„ix\a,  N)] 

~<^«) [ l2rt(jc ; a.  P>  l^)-Qn-\ix‘,oi,  p,  N)] 

where 

.  (/t+ot+/?+l)  (n+g+l)  (  N-n) 

i2n+a+p+l)i2n+ct+p+2) 

„  n  (n+P)  {n+(x+P+N+ 1 ) 

i2n+OL+P){2n+tx+p+l) 

Difference  Equation: 

-mQn(x;a,  P,  N)^Bix)[Q„ix+l;ot,  /?.  N)-Q„{x;ct,  p,  N)] 

-D(x)  [  Q^{x ;  a,  p,  N)-  Q^ix-l ;  a,  p,  N)] 

where 

fl(jf)=.(iV-jc)(a+l+x) 

Dix)  =  x{N+l+p-x) 

^.(/»)  =  n  {n+a+P+ 1) 

Rodrigues'  Formula: 
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Table  9.  JACOBI  POLYNOMIALS 


Symbol:  fJ-V) 


Interval:  [—1,1] 
Weight:  (1 -af)‘(l  + 


Standardization:  I**' '*(1)  “  ^  ^ 


Norm: 


f '  fi--'  (1  -  (1 + .y  .  j::flE(»±2±i)n2±^ 

L^«  wj  (1  u+-itr<«  (2„+a+;j+i))i!r(n+i«+/;+i) 

*'-1 

.  -  .  V  («+!)«  ^  f-w.  w+g+^-Hl  ■  1-jf  1 

Hypergeometnc  Senes:  Pr"{x)» — jjj — jf,  I  ’  2  J 

Recurrence  Relation:  2  («+!)  («+«+/!+ 1)  (2rt+a+^) 

« (2rt+a+/?+l)  [  (In+a+fi)  (2«+a+/H-2)  x+a^-p^  ] 

-2  (/»+«)  («+/?)  (2«+«+/?+2)  Ptf\x) 

DiiTerential  Equation: 

(l— +[  —  a  —  ia+p+2)  jr]y+rt  (n+a+/l+l)>' *  0 ,  y^Pji  '^\x) 
Rodrigues'  Formula: 

2%!  ??'V) «(-!)'' (i-^r(i+^r^  -^[(1 

Generating  Function: 

( 1-/+V 1  - 2„ + y  V (i+t+j I  -  it( + y )-" 

sj{l-2xt+t^)  ^ 

Explicit  Expression:  P^*'  ^(x) «  2""  S  (n-f) 


Table  W.  LAGUERRE  POLYNOMIALS 


Symbol: 

Interval:  [0,  oo) 

Weight:  x^e-*,  a>-l 

(-I)* 

Standardization:  «■  — ; —  .. 

nl 

Norm: 


Jfi  w* 


Hypergeometric  Series: 


Recurrence  Relation: 


(n+l)  I'tiW  -  [  ( 2n+«+l  )-l]  CW  -(n+a)  L”,U) 

DifTerential  Equation:  xy"  +( a+  1-.t)  y’+ny  =  0 ,  =  CV) 

Rodrigues'  Formula: 


(«)/ 


(a) 


&\x\  a _ I _ ^  r  1 


Generating  Function: 


/•■O 


/CoO 


Explicit  Expression: 


Table  11.  HERMITE  POLYNOMIALS 


Symbol:  HJix) 

Interval:  (-00,00) 

Weight: 

Standardization:  //.(I)** 

Norm: 

•'—00 

Hypergeometric  Series: 

H„ix)  -  (2jc)"  '  (j«+l)/2 ;  -  ^  j 

Recurrence  Relation: 

DifTerential  Equation: 

y"  -^Ixy  +  2ny^0 ,  y^H„{x) 

Rodrigues'  Formula: 

Generating  Function: 

H  ~ 

rt»0 

Explicit  Expression: 

//„(jr)=»rt.  ^  (-1) 

kmO  '• 
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Table  12.  GEGiyBAUER  (Xn^TRASPHERlCAL)  POLYNOMIALS 

Symbol:  CXx)  (orP.'V)).  ^>“1/2 


Interval:  [-1,1] 
Weight:  (1-x*)'-*'* 
Standardization: 


K  d^d) 


r„(x),  rt-o,  1,2,, 


Norm: 


J  LW  WJ  d  <«  (/,+i)n!r(n+2;.) 


Hypergeometric  Series: 


QW-  2^1  [  ^+1/2  •  2  J 


Recurrence  Relation: 

(«  +  1)  C<5,W  -  2  (n+;)  X  dV)  -  (n+22-l)  d/j.W 
Differential  Equation: 

(l-jc^)  j;"-(2A+l)  xy+n  {n+2X)y  -  0 ,  >>  - 

Rodrigues'  Formula:  -  1 ,  C/V) "  2^.JC, 

2"  «!  (2+1/2),  (l-x’)'-'"  C^'Vx)  -  (-1)"  (22).  [  (l-x’)"*'-'"] 


2d,V)("-(l-2xt+/’)-‘ 

naO 


Generating  Function: 


Tabh  13.  CHEBYSHEV  POLYNOMIALS  OF  THE  FIRST  KIND 


Symbol:  TJix) 

Interval:  [-1,1] 

Weight:  (l-;r>)-''» 
Standardization:  r,(l) « 1 


Norm: 


Hypergeometric  Series: 


."i"-.  V-] 


Recurrence  Relation:  "* 

Differential  Equation:  (1  -  x*)  -  .r/  +  »  0 ,  j/ »  r,(x) 

Rodrigues'  Formula; 

2-r(n+l72~^K'""^  ^ 


Generating  Function; 


2xf+r‘ 


,  -1  <x<  1 ,  I  r  I  <  1 


Explicit  Expression: 


TJx) »  cos(rt0)  with  x«cos0, 


.  W  -  f  E  ( -D*  <.^r^ 


Table  14.  CHEBYSHEV  POLYNOMIALS  OF  THE  SECOND  KIND 


Symbol:  U„{x) 

Interval:  [-1,1] 

Weight:  (1— 

Standardization:  t/,(l)  =  «  +  1 


Norm: 


Hypergeomctric  Series: 


Recurrence  Relation:  =  2xU„{x)  - 

DifTerential  Equation:  (1  -  x*)  -  3xy'  +  n  {n+2)y  =  0 ,  >^  =  U^{x) 

Rodrigues'  Formula: 

”  (l-,vY^’2'”^‘ r(rt+3/2)  ^ 


Generating  Function: 


- T— T’  IH<1 

i-ivz+G 


Explicit  Expression: 


3in(A7  +1)0  o 

- - jr  =  cos0,  0^9  <n 

"  sin0 
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Table  15.  LEGENDRE  (SPHERiCAL)  POLYNOMIALS 


Symbol:  P„{x) 

Interval:  [-1,1] 

Weight:  1 

Standardization:  P„(l)*=l 
Norm: 

Hypergeometric  Series:  P»(Jf)  “jfiF  J 

Recurrence  Relation:  («+ 1)  =  (2rt+l)  x  PJ[x)  -  n  P^iix) 

Differential  Equation:  (1  —  j/'  -  lx/  +  n  (/i+l)^'  =  0 ,  =  P,(jc) 


(-1)"  cT 


T  n\  dxT 


[(1-^)1 


OQ 

^P„(x)r''  =  (l-lxr+/^)"’/\  -l<x<l,  IH<I 

n»0 

Explicit  Expression:  P„(:c)  l’  ( -1)*  (^) 


V.  APPLICATIONS 


Orthogonal  polynomials  and  special  functions  in  general  have  been  studied  exten¬ 
sively  in  mathematics  and  other  fields  since  the  eighteenth  century.  Presented  below  are 
a  few  of  the  traditional  applications  of  selected  classes  of  orthogonal  polynomials.  Our 
first  few  applications  come  from  numerical  analysis. 

A.  ECONOMIZATION  OF  POWER  SERIES 

Economization  of  power  series  is  a  technique  used  to  reduce  the  degree  of  a 
polynomial  approximation  to  a  given  function. 

The  maximum  norm  (or  U"  —  norm)  for  a  continuous  function  on  a  compact  inter¬ 
val  Cc,i»]  is  defined  as 

||/||^=tnax  |/(jr)|. 

This  norm  is  not  induced  by  an  inner  product,  but  nevertheless  has  many  uses  in  applied 
mathematics,  including  numerical  analysis. 

The  minimax  property  of  the  Chebyshev  polynomials  states  that  of  all  n*  degree 
monic  polynomials  (i.e.,  leading  coefficient  1),  2'""r„(jf)  has  the  smallest  maximum  norm 
on  [  —1,1]  (Ref.  23;  p.  106).  The  justification  for  this  statement  is  deferred  until  Section 
C.  Hence  the  best  approximation  in  the  maximum  norm  to  the  function  xr  on  [  -1,1] 
by  a  function  of  lower  degree  is  f(x)  ^  X' —  2'-^TXx).  So,  given  a  function  and  a 
polynomial  approximation  to  that  function  (e.g.,  from  a  Taylor  series  expansion),  suc¬ 
cessively  replace  the  highest  powers  xr  with/;(jr)  to  obtain  a  polynomial  approximation 
of  lower  degree  [Ref.  23:  p.  125]. 

Example  1:  Let  f{x)  =  sin  x.  The  Maclaurin  scries  for  this  function  is 


sin  jc 


on 

-E 

*=0 


( 2*+l 

(2A-f-  1)! 


and  this  series  is  convergent  for  jc  e  R.  If  truncated  after  the  jr*  term,  the  polynomial 
approximation  for  the  function  is 


sin  x^x  — 


3! 


+ 


5! 
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with  a  maximum  error  of  0.0002  for  xe  [-1,1].  We  use  the  fifth  degree  Chebyshev 
polynomial  Ts(jc)  »  1 6jr*  —  20x’  +  Sx  to  obtain 

^  5x^ _ 5x 

5!  ~  480  1920 


with  an  error  not  exceeding 

max  I  r5(x)  I _ i_ 

1920  1920 


«  0.00052 


in  [  -1,1].  Thus  the  approximation 


smx; 


iii 

480 


75  ^3 

1920  J  1920  480 


has  an  error  whose  magnitude  in  [  -1,1]  does  not  exceed 

0.00052  +  0.0002  =  0.00072. 


Compare  this  with  the  maximum  error  of  0.00833  for  the  Maclaurin  series  which  is 
truncated  after  the  x^  term.  For  a  cubic  polynomial  approximation  of  sinx,  the  "econ¬ 
omized"  polynomial  has  a  maximum  error  that  is  significantly  smaller  (less  than  one 
tenth)  than  that  of  the  truncated  Maclaurin  series. 

The  next  three  applications  illustrate  the  usefulness  of  the  zeros  of  orthogonal 
polynomials.  The  first  two  come  from  numerical  analysis,  the  third  from  a  problem  in 
electrostatics.  We  begin  with  a  preliminary  discussion  of  a  fundamental  technique  from 
numerical  analysis. 

B.  POLYNOMIAL  INTERPOLATION 

Polynomial  interpolation  is  a  method  of  approximating  a  given  function  with  a 
polynomial  that  matches  {interpoiates  )  the  function  at  specified  points  (called  nodes  or 
abscissae)  x,, ...  ,x,  (Ref.  24:  p.  497].  Given  a  function /(x)  and  n  distinct  nodes  in  a 
compact  interval  [a, 6],  there  is  a  unique  polynomial  of  degree  (n  —  1)  that  passes 
through  the  points  (x, ,  f{x,)),  I  ^i^n. 

For  each  /  =  1,  2, ... , «,  define  a  polynomial  of  degree  (w  -  1)  by 
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Clearly,  for  each  A  «■  1, 2, 7r,^(jr*)  -0  if  k^i  and  *  1,  that  is,  ?t,^(af*)  ■=  <5,,*. 
(We  may  thus  equivalently  express  n,Jx)  =  ljix)l{(x-xi)  /,'(jc,)),  where  l„{x)  =  fl  (x-x^)  is 

kml 

the  unique  momc  polynomial  of  degree  n  that  vanishes  simply  at  each  node  x*. 
The  Lagrange  interpolating  polynomial  is  now  given  by 

n 

/a| 

Clearly,  =/(jr*)  for  /c  =  1, 2, ... ,  n,  and  uniqueness  of  the  degree  (/» —  1)  interpolat¬ 
ing  polynomial  is  guaranteed  by  the  Fundamental  Theorem  of  Algebra. 

C.  OPTIMAL  NODES 

In  this  section  we  address  the  issue  of  estimating  the  maximum  size  of  the  interpo¬ 
lating  error  ||/-L'IL.  Assuming /is  suitably  differentiable  in  [a,i5>],  it  can  be  shown 
(Ref.  25:  p.  188]  that  there  exists  a  value  K,  (which  depends  on/*” in  such  that  for 

any  x  in  [a,f>]. 


H 

where  l^x)  “  FI  (jf  “  ■*■()•  (Note  again  that  l„ix)  is  a  monic  polynomial  of  degree  n  which 
1-1 

vanishes  at  the  nodes.)  This  implies  that  ||/—  L;  |L  is  minimized  by  making  the  optimal 
choice  of  nodes  jc„  ... ,  x,  in  [a,6]  which  minimizes  ||  /,  |L.  Surprisingly  perhaps,  this  hap¬ 
pens  precisely  at  the  zeros  of  the  Chebyshev  polynomials  T„{x),  scaled  to  the  interval 
[a.f)].  We  sketch  the  reasons  below. 

For  simplicity,  we  take  our  interval  of  interest  to  be  [  -1,1]  instead  of  without 
loss  of  generality.  The  transformation 


or  equivalently, 


is  a  one-to-one  continuous  mapping  between  the  intervals  [  — 1,1]  and  [a,i]. 

Now  let  X,  be  such  that  TJix,)  =  0,  /  =  1, 2, ... ,  n  (sec  Chapter  I,  Section  A.3),  i.e.,  let 
Ffona  the  recurrence  relation  (Equation  (4)  in  Chapter  I),  TJ[x)  has  a 
leading  coefiicient  of  2"“'  ;  hence  this  /,(x)  is  monic.  Moreover,  from  the  definition 
TJi  cos  6)  -  cos  nB,  it  follows  that 


(1) 

where  y,  ,  /  =  0, 1, ... ,  n.  Now  suppose  /rt,(jf)  is  another  monic  polynomial 

of  degree  n  such  that 

(2)  = 

Combining  (1)  and  (2)  we  see  that  we  must  have  for  i  *  0, 1 . n, 

J'/)  <  Uyd  =  “Tib"  ^ 

(3) 

>  Iniy^  =  ^;;zr  odd. 

Thus,  the  polynomial  p„{x)  =»  n\(x)  -  l„ix)  has  degree  at  most  («  -  1)  (since  it  is  a  difier- 
ence  of  two  monic  polynomials  of  degree  n)  w;ith  at  least  n  zeros,  one  in  each  interval 
(y  .J'z+i).  0  ^  ^  rt  —  l,  by  (3).  This  contradicts  the  Fundamental  Theorem  of  Algebra, 
and  so  no  such  polynomial  m„(x)  satisfying  (2)  exists.’  Hence  the  choice  of  nodes  or,  de¬ 
termined  by  /„(x)  =s  -z-T-  T„(x)  minimizes  ||  /„  ||„  =»  max  |  /„(x)  |  over  all  possible  monic 
polynomials  /„(x),  and  therefore  over  all  possible  choices  of  interpolating  nodes  jc,. 

It  should  be  emphasized  that  this  Chebyshev  interpolation  allows  an  a  priori  error 
bound  for  all  x,  but  is  not  always  best  possible  for  every  jc  using  other  interpolation 
schemes.  (For  example,  if  equally  spaced  nodes  jr„ ... ,  jr,  arc  used,  then  trivially,  the  er¬ 
ror  fix,)  —  mx)  -  0,  even  if  ,y  =  x,  is  not  a  Chebyshev  zero.)  Chebyshev  zeros  arc  opti¬ 
mal  when  one  has  freedom  in  the  choice  of  nodes.  For  details  regarding  the  practicality 
of  Chebyshev  interpolation  and  some  asymptotic  results,  sec  (Ref.  24]. 
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D.  GAUSSIAN  QUADRATURE 

Quadrature  formulas  are  used  in  that  area  of  numerical  analysis  concerned  with  the 
approximate  integration  of  a  function /(jc)  against  a  weight  function  w{x)  >  0  on  an  in¬ 
terval  (a,b),  when  the  explicit  evaluation  is  intractable. 

An  interpolatory  quadrature  is  such  a  rule  that  uses  interpolating  polynomials,  such 
as  Lagrange  polynomials: 


J  fix)  w(jf)  «  J  liix)  w(jf)  dx, 

where  L{{x)  =  S/W  xfx),  as  described  above.  This  can  be  rewritten  as 
1-1 

(1)  I  /(jf) 

*'«  i-i 

where  the  weights  w,  are  given  by 

(2)  w,  =  f  w(jf)  dx. 

•'ll 

Thus,  interpolatory  quadrature  is  basically  a  weighted  sum  of  the  function  values /(xj 
at  the  nodes  jr, ,  /  =  1, 2, ... ,  n,  as  are  numerical  integration  recipes  such  as  Simpson's 
Rule  and  the  Trapezoidal  Rule. 

For  specified  nodes  Xi, ... ,  jc,  e  [a,6],  tlie  n  weights  w,, ... ,  w,  computed  in  (2)  for  the 
quadrature  (1)  will  be  exact  for  polynomials  /e  F,_i[a,6],  but  we  can  do  better.  [Ref. 
25:  p.  2361 

In  Gaussian  quadrature,  we  ask  for  the  location  of  the  n  nodes  jc„  ... ,  jc,  as  well  as 
the  n  weights  w,, ... ,  w,  in  order  for  the  quadrature  rule  (1)  to  be  exact  for  polynomials 
/e  Pin-Sfiyb'].  At  first  glance,  this  seems  to  be  an  extremely  complicated  computational 
problem,  but  the  solution  falls  out  simply  when  the  theory  of  orthogonal  polynomials 
is  applied. 

Let  { p„{x)  }“o  be  the  class  of  polynomials  orthogonal  with  respect  to  the  weight 
function  w(jf)  on  say,  by  Gram-Schmidt,  and  let  L[[x)  be  the  Lagrange  polynomial 
that  interpolates /e  at  the  zeros  of  p„{x),  so  that/(.r)  —  L[{x)  e  van¬ 

ishes  at  the  zeros  also.  Since  the  zeros  of  p„(x)  arc  real,  simple,  and  lie  in  (a,b),  we  have 
the  property  that 
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fix)  -  l{(x)«p„(x)  <f^i(x). 


where  e  iViCfl.A].  Hence,  by  orthogonality, 

[  /(x)  M<af)  </x-  f  l{(x)  H<jf)  </jr  -  [  p^{x)  »v(jr)  dx  -  0. 

•^a  '^a  •'a 

We  have  thus  shown  that  the  resulting  n-point  Gaussian  quadrature  rule  is  exact  for 
f  ^  Pbt-\La,b2.  Moreover,  by  a  theorem  of  Sticltjes,  IT  fix)  is  continuous  on  a  finite  inter¬ 
val  [o,i],  then 


.‘is,  f 


dx. 


By  applying  the  Christofiel-Darboux  formula  (Equation  (9)  of  Chapter  III  with 
notation  from  Chapter  III,  Section  B)  and  using  (1),  it  is  possible  to  derive  an  alternate 
expression  for  the  weights  (2): 


KPfixdPn-^liXi) 

which  are  referred  to  as  the  Christq,Jel  numbers.  (Ref.  8  ] 

Thus  GausS'Jacobi,  Gauss-Chebyshev,  and  Gauss- Legendre  are  the  names  given  to 
Gaussian  quadratures  involving  the  weight  functions  and  orthogonal  polynomials  from 
the  Jacobi,  Chebyshev,  and  Legendre  classes,  respectively. 

E.  ELECTROSTATICS 

The  zeros  of  the  Jacobi  polynomials  play  an  interesting  role  in  a  problem  of  Sticltjes 
concerning  electrostatic  equilibrium.  In  this  problem,  fix  "masses’  of  positive  charge  a 
and  at  the  points  x  =*  1  and  jf  =>  - 1 ,  respectively.  Then  place  n  point  masses  of  positive 
unit  charge  in  the  interval  ( -1,1)  so  that  they  are  free  to  move.  These  interior  masses 
are  now  subject  to  a  "logarithmic  potential",  that  is,  a  repelling  force  that  is  proportional 
to  the  logarithm  of  the  distance  separating  them.  The  problem  is  to  detcrnainc  the  dis¬ 
tribution  of  the  point  masses  jf, ,  /  =  1, 2, ...  ,/i,  when  the  system  is  in  equilibrium. 
Mathematically,  this  is  equivalent  to  maximizing  the  force  function 
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yi  n 

<3)  ... ,  jr„)  -  o  log(l-x,)  +  /?  2]  log(i+Jf/)  +  ^  log  1  ^i-Xj  | . 

i—l  /—I 

The  logarithmic  terms  give  the  restrictions  jr, 9^“— 1,  1,  and  x^^Xj  for  ii>j.  This 

otherwise  continuous  function  gives  the  equilibrium  points  by  setting 


(4) 


dF 


Sxi 


-0,  /=  1,2, ... ,rt 


and  solving  for  jr, . jr,.  To  solve  this  system  of  n  nonlinear  equations  in  n  unknowns, 

Stieltjes  introduced  the  polynomial 


and  reduced  (4)  to 


1  » 

2  Pn'W  ^  l+^t, 


a 

1-jr, 


=  0,  1-1,2, 


n 


which  becomes 

(l-ar'>/(.r,)  +  [  2/?  -  2a  -  (2a  +  2/1)  .t,  ]  -  0 

for  i »  1, 2, ... ,  n. 

Since  the  polynomial 

( 1  -  x^)  p„%x)  4-  [  2/?  -  2a  -  (2a  +  2/1)  x  ]  pM 

is  of  degree  at  most  n  and  vanishes  at  jf  =  jf„  i  —  1,  2, ... ,  n,  it  can  be  set  equal  to  a  scalar 
multiple  ^pXx)  which  also  vanishes  at  these  points.  Hence 

(1  -  .r^)  p„“{x)  +  [  2)?  -  2a  -  (2a  +  20)  x  ]  p„'ix)  -  Ap„ix)  =  0. 

Attempting  a  power  series  solution  to  this  second-order  dilferential  equation  leads  to  the 
observation  that  polynomial  (i.e.,  terminating)  solutions  exist  if  and  only  if 

X  =  —n[n  +  2a  +  2/1  —  1). 
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Rearranging  terms  yields 


(1  -  x^)  p„%t)  +  [  (2/1  -  1)  -  (2b  -  1)  -  ((2b  -  1)  -H  (2/1  -  1)  +  2)  ;r  ]  pM 
+«(«  +  (2a  -  1)  +  (2/S  -  1)  +  I)  p„{x)  -  0 


which  is  the  differential  equation  for  Thus  the  equilibrium  positions  of  the 

unit  charges  occur  at  the  zeros  of  this  Jacobi  polynomial.  [Refs.  8,  13:  p.  140] 

The  zeros  of  the  Laguerre  and  Hermite  polynomials  can  be  developed  as  solutions 
to  similar  electrostatic  equilibrium  problems. 


F.  SPHERICAL  HARMONICS 

We  investigate  another  application  of  orthogonal  polynomials  to  problenis  in 
mathematical  physics.  In  Cartesian  coordinates  (jf,.v,  z),  the  Laptacian  operator  of  a 
function  u(x,y,2)  is  defined  as 


.  d^u 

dx^  6/  dz^  ‘ 


When  converted  to  spherical  coordinates  (r,  9, 0),  this  operator  acting  on  a  function 
u(r,  6, 0)  becomes 


(5) 


V^u 


( sin  9  0^)0  + 


1 

sir?9 


U0$ 


where  subscripting  with  r,  9,  or0  denotes  partial  differentiation  with  respect  to  that 
variable.  A  function  u  is  said  to  be  harmonic  in  a  region  D  if  in  that  region  it  satisfies 
Laplace's  equation: 

*  0. 


In  particular,  if  the  boundary  of  D  is  the  unit  sphere  centered  at  the  origin,  then  using 
the  method  of  separation  of  variables,  the  Legendre  polynomials  arise  naturally  as  part 
of  the  solution. 

In  separation  of  variables,  we  assume  that  the  solution  will  be  of  the  form 
(6)  u{r,  9,  <f}) -f{r)  g{0) 

Substituting  (6)  into  (5)  and  dividing  through  by  u,  we  obtain 
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(7) 


/  g  sinO  ^  A^^-O. 


The  first  term,  depending  only  on  r,  must  reduce  to  a  constant  which  we  write  as 
v(v+l).  Substituting  v(v+l)  into  (7)  and  multiplying  through  by  sin*0,  we  have 

v(v+l)  sin^a  +  (sin0^#)j+  ^ 


Now  we  see  that  the  third  term,  depending  only  on  0,  must  also  reduce  to  a  constant; 
call  it  -m*.  Substituting  -/«*  for  this  term  and  simplifying,  we  obtain 

(8)  sin^0g^^+  sind  cos0g^  +  [v(v+l)sin*0  — 

By  the  change  of  variable  x  »  cos  6,  (8)  becomes  via  the  Chain  Rule 

For  problems  that  are  radially  symmetric,  the  0  dependence  can  be  removed  by  setting 
m  «  0,  leaving 

(9)  (l-x’)-i-2j:-^+v(v+l)*-0. 

ax 


When  v  is  a  positive  integer  n,  (9)  is  recognized  as  the  differential  equation  for  the 
Legendre  polynomials,  [Ref.  26:  pp.  210-213] 

An  alternate  approach  to  this  problem  (by  the  method  of  images)  uses  Greens 
functions.  As  motivation,  consider  the  Dirichlet  problem  for  the  unit  circle  in  the  plane, 
which  involves  finding  a  harmonic  function  u(r,  Q)  in  the  unit  disk  that  takes  on  pre¬ 
scribed  function  values  f[d)  on  the  boundary  r  =  1.  The  solution  is  given  by 


where 


I-r' 


I— 2r  cos  y+r^ 
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is  the  so-called  Poisson  kemtl  fbr  this  probkot 

Similarly,  the  solution  to  the  IMckkt  ike  mmi  tfkeee  m  ft* 

on  r«  I 


can  be  expressed  as 


where 

(10) 


PM^ 


A=±. 


[  l~2rcosr+r*]*'* 


and 


cos  y  «  cos  9  cos  9'  -f-  sin  9  sin  9*  cos(^  •> 

Writing  ( 10)  in  terms  of  simpler  functions 

I 

y/l'-2r  cos  y+r* 

we  note  the  appearance  of  the  generating  function  for  the  Legendre  pol>'nomials,  with 
X  -  cos  y.  Therefore  the  Legendre  polynomials  are  again  part  of  the  solution.  [Ref.  27: 
pp.  87-89] 

Laplace's  equation  can  also  be  solved  in  a  higher  dimensional  setting.  Let 
X  =  (jf„  ... ,  jf,)  6  R'.  Note  then  that  we  may  .  wite  x  -  r\,  where 
r  “  II X II  *  ^/xf  +  ...  -H  xj  and  ^  i, ... ,  ^,)  e  R'  is  a  unit  vector.  A  polynomial 

hJix)  is  said  to  be  homogeneous  of  degree  n  if  A,(^.x)  =  A"  /i,(x). 

The  Laplacian  operator  acting  on  a  function  u(x)  in  p  dimensions  is  defined  as 


PM  ■■ 


I 


cos  y+r 
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As  before,  a  harmonic  Ainction  is  one  which  satisfies  Laplace's  equation  V|u  >■  0.  We 
now  seek  homogeneous  harmonic  polynomials  of  degree  n  in  x  e  R'.  It  can  be  shown 
that  there  are  exactly 

-3 
1 

linearly  independent  such  solutions,  and  they  can  be  characterized  by  Gegenbaucr  (or 
ultraspherical)  polynomials  A  «  1, 2, ... ,  N.  The  general  solution  is  given  by 

Mx)  -  hM)  -  -SnCO. 


2n  +  p  —  2 


n 


(■" 


where  the  spherical  harmonic 

km\ 

with  ii  suitably  chosen  unit  vectors.  Note  that  if  p  *  3,  then  these  reduce  to  the 
Legendre  polynomials  found  earlier.  (Ref.  6:  pp.  168>183] 

G.  GENETICS  MODELING 

Karlin  and  McGregor  gave  an  interesting  application  of  the  dual  Mahn  polynomials 
to  a  model  in  genetics.  In  this  continuous  time  Markoff  chain  model,  the  dual  Hahn 
polynomials  R,(  ^.(x) ;  y,  d ,  iV )  arise  in  the  transition  probability  function  for  tlic  process. 

The  setting  for  the  model  assumes  N  gametes  of  type  a  or  A  and  gives  a  random 
fertilization  scheme.  The  population  of  either  type  of  gamete  is  affected  by  both  the 
fertilization  process  and  a  mutation  process  whereby  a  gamete  resulting  from  a  mating 
can  mutate  into  the  other  type. 

By  considering  the  conditional  probabilities  for  an  increase  in  population  size  of 
both  gamete  types,  a  stochastic  process  is  defined  which  is  a  classical  birth  and  death 
process.  The  transition  probability  function  for  this  last  process  is  then  cast  in  terms 
of  the  dual  Mahn  polynomials.  The  interested  reader  is  referred  to  (Ref.  14]  for  the  de¬ 
tails. 
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VI.  BASIC  EXTENSIONS 


A.  BASIC  HYPERGEOMETRIC  SERIES 

In  this  chapter,  we  extend  the  structure  and  some  of  the  results  in  Chapter  IV  to  a 
more  general  level.  This  extension  is  accomplished  by  introducing  a  new  parameter 
called  the  base  to  the  hypergeometric  series.  The  base  q  was  used  by  Heine  in  a  series 

.  (l-v»)(lV)  (!-,«)  (1-0  (lV)(l-,H.)  ^ 

where  c  0,  —1,  -2, ...  (Refs.  28,29].  This  series  converges  absolutely  for  I  x  |  <  1  when 
I  9  I  <  1  by  the  Ratio  test.  Since 

(2)  lim  \  »  a, 

v-*i  1-? 

we  sec  that  the  series  in  (1)  tends  termwise  to  the  ordinary  hypergcometric  scries  as 
^  1.  Thus  Heine's  series  is  called  the  basic  hypergeomeiric  series  or  the  q- 
hypergeometric  series,  [Ref.  12;  p.  3] 

The  q'Sh{fied  factorial  is  the  basic  extension  of  the  shifted  factorial  introduced  in 
Chapter  II  and  is  defined 

fl,  rt-0 

(3)  (« :  q)n  -  »  1,  2 . 

The  ordinary  shifted  factorial  is  recovered  by  applying  (2)  and  (3)  in  the  limit 

(?  •  4)n  ,  X 

v-»i  (l-^) 

We  also  define 


(« :  =  n 
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a  form  we  will  see  in  later  results.  Since  the  infinite  produa  diverges  when  botli 
and  I  q\  ^  1,  we  will  assume  ]  <7  |  <  1  whenever  {a ;  q)^  appears  unless  otherwise  stated. 
[Ref.  12:  p.  31 

Generalizing  the  basic  hypergeometric  series  above,  we  define  the  ,0,  basic 
hypergeometric  series  (or  ,<t>,  series) 


jA  ,  ..1  ^  (^t  I  4')rt  —  (^r  I  ?)n  ff  IN'*  ,0  1" 

where  (”)"«(«— 1)/2  and  ^#0  when  r>j+l.  We  require  that  the  parariieters 
bi, ... ,  b,  be  such  that  the  denominator  factors  in  each  term  of  the  series  are  nonzero. 
Since 


(<7  =  0 ,  /I  =  m  +  1,  m  +  2, ... , 

a  ,(t>,  series  terminates  if  one  or  more  of  the  numerator  parameters  is  of  the  form  </■"  for 
m  *  0, 1, 2, ...  and  q^O.  When  r  =  s  +  1,  the  expression  in  (5)  simplifies  to 

ra„...,a,+,  1  y  (fli ;  q)„ ...  ;  q)„  „ 

”'^’L  *1 . *,  J“„4s  («;?). (*1 ; 9). 

Note  that  in  a  basic  hypergeometric  series  c„  x",  the  ratio  c„+,/c,  is  a  rational  function 
of  q”.  [Ref.  12:  p.  4] 

Using  the  q  —  shifted  factorial  (3)  and  (4),  we  can  define  basic  extensions  for  many 
of  the  functions  and  formulas  introduced  in  earlier  chapters.  We  note  that  often  there 
is  more  than  one  way  to  extend  a  result;  examples  will  be  given  shortly. 

The  q-gamma  function  is  defined  by 

(6)  r,W  =  (l-?)'"' .  0  <  q  <  1. 

(<7  i  ^)oo 

Gosper  showed  that 

(7)  Urn  r,(x)  =  r(.T) 

for  0<<ir<  1  [Ref.  30:  p.  109].  The  structure  of  the  gamma  function  extends  as  well. 
For  instance,  the  formula 
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can  be  reduced  using  (6)  to 


r(;f+l)-xr(jc). 


With  (6)  we  can  deflne  the  q~beta  function 


r,U)r,(>>) 

r^(x+y) 


which  by  (7)  tends  to  Bix,y)  as  ^  1-. 

The  q-binomial  coefficient  is  defined  for  integers  n  and  k  by 


r«"|  (? ;  9)n 


where  A  »  0, 1, ... , «.  For  nonintegral  a  and  /S,  we  define 

Jj<i  iq ;  q)oo  (7”^’ ;  9)00 

r^(a+l) 

"  r,(/i+i)r,(a-/?+i)  • 


The  q'-binomial  theorem  is  then 


/c-0 


where  /?  =  0, 1, 2, ... .  [Ref.  12:  p.  20] 

The  next  two  expressions  arc  basic  extensions  of  the  Chu-Vandermondc  formula 
(Chapter  III,  Section  F.l)  and  are  both  known  as  the  q-Chu-Vandermonde  formula: 


i^\ 


q~’',  b  c<7"  1  {cib ;  q)„ 

C  ■’'  *  J“  (c:?)„ 
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[Ref.  12:  p.  11].  These  fonns  can  be  shown  to  be  equivalent  by  reversing  the  order  of 
summation. 

Jackson  introduced  the  general  form  of  the  q-integral 

J  Vw  -  JVw  "V  -  [Vw  rf,' 


where 


OO 

f/W-i,'-- <<!-?) 

The  ^  —  integral  defines  the  measure  which  is  a  natural  object  for 
q  —  defined  functions  (Ref.  12:  p.  19).  All  of  the  functions  and  formulas  developed  above 
play  an  important  role  in  generalizing  the  ordinary  orthogonal  polynomial  classes. 

B.  BASIC  EXTENSIONS  OF  ORTHOGONAL  POLYNOMIALS 

In  this  section,  we  present  the  q  —  analogue(s)  of  selected  classes  from  earlier  chap¬ 
ters.  By  using  formulas  such  as  those  presented  in  the  previous  section  together  with 
methods  based  on  those  outlined  in  the  preceding  chapters,  it  is  possible  to  derive  the 
recurrence  relations,  difference  equations,  and  Rodrigues'  formulas  as  well  as  many 
other  identities  satisfied  by  these  q  -versions.  As  mentioned  in  the  previous  section,  the 
q  —  extension  of  a  function  is  not  necessarily  unique;  however  the  last  two  classes  listed 
are  especially  portant. 

1.  Continuous  q  — Hermite  Polynomials 
a.  Definition 


(<y 9)n 

(<7 ;  *i)k  (<7 ;  ^)n-k 


^Kn~2k)9 


where  jc  =»  cos  6, 


b.  Orthogonality  Relation 


where  x  >■  cos  d,  [Ref.  12:  p.  188] 

2.  Discrete  q  — Hennite  Polynomials 
a.  Definition 


[»/2] 


Y,  71^ 

“•  I  /I*  •  /? 


;  <7)n 


Orthogonality  Relation 


r 

•'-1 


//^(x ;  »  q 


GK:,) 


n 


where  \liix)  is  a  step  function  with  jumps 


I  jt  I  (jtV  w  LUw ). 
2 


at  the  points  x  =«  ±  ^ ,  y  «  0, 1, 2, ... .  [Ref.  12:  p.  193] 
3.  q  -Laguerre  Polynomials 
a.  Definition 
For  a  >  —1, 


/(«)/  \  '» k\^  .  /»+«+i"l 


Continuous  Orthogonality 


f  ”  i}^\x .  - -  r(a+l)r(  at)  (<7“-^' ;  <7)^ 

Jo  r,(-a)  U ;  (?)„  <7 
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c.  Discrete  OrtkegomUity 


2 


(-c(l-^)/;i?)^ 


where 


;  ^)oo  (-(i-?)/(c?*) ;  yL 

;  ^)oo  (“c(J“^) ;  9)^  (“  T  i  9)00 


[Ref.  12:  pp.  194-195] 

4.  Little  f  -  Jacobi  Polyaomiab 
a.  Dtfinitien 


p,(jf  \a,b',q) 


'2^1  ^ 


abq 

aq 


,n+l 


W.  xq 


b.  OrthegeiiaBty  Rdmtien 


^  Pm{9*  ;«.*;?)  pn{9^ ; ;  ?)  ^^9) 

i9 ;  ?)« (l-oA^)  ibq ;  <?)«  ;  q)^ 


iohq ;  (l-ab^^"'*'*)  (aq ;  q)„  {aq  \  <7), 


m.n 


where  0  <  ^ ,  aq<l  [Ref  12:  p.  166] 
5.  Big  q  -  Jacobi  Polynomials 
a.  Definition 


n  /  t  \  ,  9  .<*»?.  Jf 

P;,(jr ;  a,  b,  c ;  -  3<^2  ;9,  17 

aq,  cq 
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A.  OrthpgonaUty  Rehtion 


r 


P„ix ;  a,b,c  ;  q)  P^ix  ',a,b,ciq) 


ixlaiq)^ixlciq)^ 

ix\q)^ibxlc\q)^ 


i<^bq ;  q\  (l-aA^*  +‘)  {aq ;  q\  {cq ; 


where 


, ,  r*’  {xfa\q)^{xic\q)^ 

Jc,  (x‘,q)^ibxlc;q)^ 

aq  (t-<7)  (g ;  y)oo  (c/a ;  y)oo  (q?/c ;  ;  q)^ 

(oq ;  q)^  (bq ;  q)^  {cq ;  q)^  {abqjc ;  q)^ 


[Ref.  12:  pp.  167-168] 

6.  q  "Krawtchouk  Polynomials 
a.  Definition 


K„{x ;  a,  N  \  q)  -  3^2  *^“^0  ^  ^  ] 


A.  Orthogonality  Relation 


ff 

v-i 


V  w,"  i «.  AT ;  f)  *;(»-' : Af :  ?)  (-a/ 

JCO 


[Ref.  12:  p.  185] 

7.  q  -Hahn  Polynomials 
a.  Definition 


m/\ 


C„(jf ;  a.  A,  yV ;  <jf)  =-  3^2  ;  <7*  '7 

L  9  J 
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b.  OrthogtnaUty  Relation 


ff 


2  CJx : iv  i ,)  a(x ; t.  JV : ,)  («»)-' 

(y ;  (l~afe^)  jbg ;  ;  q)„ 

iabq ;  (l-aA?*'"^')  (a^ ;  q)^  (q’'^' ;  q\ 


x-O 


(9 ;  4)s 

X  (-a^)"  ^(a)“^  5^^ ,  m,  n  -  0, 1, ... ,  AT 


[Ref.  12:  p.  165J 

8.  Dual  f  -Hahn  Polynomials 
a.  Definition 


Rni  M(jf) ;  A,  c,  iV ;  ^)  -  302  ’  ’  Zv  * .  1 ^ 

q  ,  bcq 


where  mW  — 

A.  Orthogonality  Relation 


N 

^  R„(M(jf) ;  A,  C,  iV ;  q)  R„in{x)  \b,c,N',q) 


XmO 


(c?~^ ;  ?)x  (^c? ;  ?)x  ;  y)x 

;  ?)x  ;  9)x  (’ ;  9)x  ^ 


ibq\qh  ibcq\q)„{q~'  ',q), 


i-bcq)-^ 

(cq~^''  6^^  ,  m,  «  -  0,  1, ... ,  iV 


(Ref  12:  p.  1661 

9.  q  -Racah  Polynomials 
a.  Definition 


(8) 


PnMx)',a,  b,  c,«/;<7)-403 


—If 


abq"'^\  q~*,  cdq 
aq,  bdq,  cq 


x+l 


where  fx(x)  —  q-‘  +  cdq"*' 
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b.  0rtb9gi$aBty  Relation 

When  aq^  bdq^  ox  cq^^q-"  for  a  positive  integer  N,  then 


(9)  2]  \a,b,c,el\  q)  pMx)  \a,b,c,d\  q)  h<x  ;  a,  6,  c,  rf ;  6„^ 

XmC 

where 


w{x',a,b,  c,d\q) 

(10)  {cdq ;  q)^  (aq ;  q)^  (bdq ;  q)^  (cq ;  y);, _ 1  -  cdq^' 

(9 1  <r)x  (o~^cdq  ;  q)^^  (b^'cq ;  (dq ;  q)j,  (abqf  *  “ 

and 


(? ;  <r)n  (l-abq)  (bq ;  q)„  (ad~'q ;  (a^c~V ;  q)„  {cdqf 
(abq ;  q)^  {aq ;  q)„  {bdq ;  q)^  {eq ;  q)^ 

jedq^ ;  ?)oo  q)^  {b^^ ; 

{a^^cdq ;  q)^  {b'^cq ;  q)^  {dq ;  q)^  (a"'6"  V*  1 9)00 


[Ref.  20:  p.  1014).  When  aq ,  bdq ,  or  cq  is  equal  to  q-’',  the  infinite  products  in  (11) 
reduce  to  finite  products.  Hence  the  orthogonality  relation  (9)  is  valid  for  all  q  provided 
no  zeros  are  introduced  into  denominator  terms  [Ref.  18:  p.  4]. 

10.  Askey- Wilson  Polynomials 
a.  Definition 


(12) 


p„(x ;  a,  d,  c,  </  I  ?)  -  a"'*  {ab  \  q)^  {ac ;  q)„  {ad ;  q)„ 

\  abcdq'^\ 
ab,  ac,  ad 


X4^3 


[■ 


ae  ,  ae 


where  x  »  cos  0  [Ref.  18:  p.  3] 

b.  Orthogonality  Relation 
For  —1  <  a,  b,  c,  d,  q<  1, 


1  r  ^ 

■57  J  ^Pm(x\a,b,c,d\  q)pj^x\a,b,c,d\  q) 


w{x ;  a,  b,c,d\  q)  dx 

yj\-x^ 
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where 


n  (I -2(2*’ -!)<(*  +  ?“) 

'*'(*  t  a,  c,  rf  I  “  /i(x,a)  A(x,&)  /i(x,c)  A(x^ 


with 


»0 


M) 


and 


(a6crf^^" ;  q)^  (abedq'^^ ;  ?)„  (y"^* ;  q)li  (abq'' ;  y)"* 

“  (ae^" ;  4)^  {adq’' ;  (Ac?" ;  q)^  (Arf?" ;  ?)^  (erf?" ;  ?)^  ‘ 


[Ref.  18;  pp.  1M4] 

We  can  establish  a  formal  connection  between  these  last  two  classes.  If  the 
?  -Racah  polynomials  are  written 


403 


a'A'?"’^‘,  ?”',  c’ttq 
a'q,  Vdtq,  c*q 


--Jt 


Jt+i 


W.  9 


and  the  following  parameter  changes  arc  made: 


(or  likewise  using  any  permutation  of  {A,c,rf}  assuming  at  least  one  of  these  is  non 
zero),  as  well  as  the  change  of  variable  ?~*  -■  ae*,  then  we  obtain 


403 


_n— I 


abedq  '  at 
ab,  ac,  ad 


» 


ae 


9 


with 


Ai(x)-?-^  +  ?^-"’c'rf' 

—  +ae~^  ■>  2a  cos  0 


97 


That  is,  the  q  -Racah  and  Askey* Wilson  polynomials  are  virtually  the  same,  difTcring 
only  in  their  parameters,  normalization,  and  variable. 

Because  the  q  -Racah  and  Askey-Wilson  polynomials  are  essentially  the 
same,  the  names  are  oRen  used  interchangeably  in  the  literature.  We  have  presented 
these  polynomials  as  distinct  classes  in  order  to  emphasize  the  continuous  and  discrete 
natures  as  expressed  in  the  orthogonality  relations. 

As  stated  at  the  beginning  of  Chapter  IV,  the  ''classical  orthogonal 
polynomials"  are  defined  to  be  those  which  are  special  or  limiting  cases  of  the  Askey* 
Wilson  (21)  or  f -Racah  polynomials  (18)  (Ref.  31:  p.  57].  We  can  now  make  this 
statement  a  hit  more  precise.  Letting 

a-f*  ,  6-^  ,  c^q^  ,  d»q* 

in  (9)  and  (10)  and  taking  the  linnit  as  7-»l,  the  ordinary  Racah  polynomials 
R,(  >l(x) ;  a, /?,  y,  5)  and  their  weight  function  as  defmed  in  Formulas  (30)  and  (32), 
Chapter  IV,  Section  C.3  are  retrieved.  A  similar  limiting  process  will  recover  any  ordi¬ 
nary  orthogonal  polynomial  class  from  its  q  -extension.  Moreover,  any  q  -extension 
previously  discussed  is  a  special  or  limiting  case  of  the  q  -Racah.  For  example,  letting 
eqmq-ff  and  d^O  produces  the  ^-Hahn,  etc.  In  this  way,  we  sec  that  the  self-dual 
q  -Racah  polynomials  encompass  all  the  previous  classes.  They  also  satisfy  three-term 
recurrence  relations,  second  order  dilTcrence  equations,  and  Rodrigues'  formulas  with 
respect  to  q  -divided  difTcrence  operators.  The  interested  reader  is  referred  to  (Refs.  18, 
12,  20]  for  details  about  these  very  rich  classes. 

C.  CONCLUDING  REMARKS 

In  Chapter  V,  we  presented  a  few  of  tlie  traditional  applications  in  which  orthogonal 
polynomials  have  played  an  important  part.  More  recently,  an  enriching  interplay  has 
developed  between  the  theory  of  orthogonal  polynomials  and  other  mathematical  and 
mathematically-related  areas. 

Efllcient  computational  methods  have  been  devised  for  determining  the  many  useful 
quantities  associated  with  orthogonal  polynomials  (such  as  their  zeros,  recurrence  coef¬ 
ficients,  etc.)  [Ref.  32:  pp.  181-216).  Various  classes  of  orthogonal  polynomials  have 
also  played  a  role  in  digital  signal  processing  (Ref.  32:  pp.  1 15-133),  quantum  mechanics 
[Ref.  32:  pp.  217-228),  and  birth,  death  processes  (Ref.  32:  pp.  229-255).  .Advances  in  the 
field  of  combinatorics  and  graph  theory  have  allowed  new  geometric  interpretations  of 
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orthogonal  polynomial  identities,  sonw  of  which  have  very  important  consequences  for 
"association  schemes"  and  the  designs  of  codes  {Refs.  33,  16, 34,  32  :  pp.  25-33,  35]. 

Physicists  have  introduced  various  versions  of  "diagrammatic  methods":  ways  of 
understanding  orthogonal  polynomials  through  3-]  and  6-j  symbols,  their  generaliza¬ 
tions,  and  accompanying  identities  by  formally  associating  them  with  pictorial  schemat¬ 
ics  representing  forces  of  physical  systems  that  conserve  angular  momentum  [Ref.  36). 
PowerAd  new  techniques  involving  "quantum  groups"  have  been  used  to  generate  new 
identities  for  some  classes  (Ref.  32:  pp.  257-292].  Further  investigation  into  the 
electrostatics  problem  discussed  in  the  text  has  led  to  the  formation  of  the  famed 
"Selberg  beta  integral"  and  its  generalizations.  It  has  yielded  to  analysis  via  the  study 
of  the  "root  systems"  of  Lie  algebras,  and  has  found  applications  ranging  from  statistical 
mechanics  to  computer  algorithm  complexity  (Refs.  30:  pp.  48*52,  32  :  pp.  311-318,  37]. 

Finally,  research  into  the  general  structure  of  q  -series  has  led  to  many  surprising 
connections,  and  is  intimately  related  to  the  many  remarkable  and  po\>*errul  number- 
theoretic  formulas  discovered  by  S.  Ramanujan,  the  famed  Indian  mathematical  prodigy 
(Refs.  30:  pp.  87-93,  38:  pp.  55-66,  39]. 

This  utility  and  promise  of  future  applications  provide  ample  motivation  and  justi¬ 
fication  for  continued  study  of  the  intrinsic  structure  of  orthogonal  polynomials. 
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